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1 • INTRODUCTION 

Differentially  encoded  phase-shift  keyed  (DPSK)  binary  signalling 
is  a commonly  used  technique  in  digital  data  transmission  systems.  There 
are,  of  course,  other  modulation  techniques  that  are  more  efficient  than 
the  DPSK  systems  in  terms  of  power  and  spectral  occupancy;  however,  the 
circuit  simplicity  and  the  accompanying  cost-effectiveness  of  the  DPSK 
system  often  becomes  the  over-riding  reason  for  its  application  in  a 
data  communication  system. 

The  choice  of  DPSK  modulation  is  particularly  appropriate  when  a 
small,  low-cost  modulator  is  needed.  A typical  application  might  be  a 
buoy  which  has  limited  space  and  limited  prime  power  for  the  conmuni cations 
system.  An  expendable  platform  such  as  a buoy  must,  of  course,  be  low- 
cost.  The  data  transmission  from  the  buoy  may  involve  a hard-limiting 
satellite  channel  for  the  purpose  of  processing  the  data  at  a shore 
facility. 

Surprising  as  it  may  sound,  however,  a complete  understanding  of  the 

error  behavior  of  DPSK  system  has  been  lacking.  Much  less  understood  is 

the  error  performance  mechanism  of  a DPSK  system  over  the  hard-limiting 

satellite  channel.  The  error  behavior  of  DPSK  system  over  a "linear 

channel"  has  always  been  assessed  on  the  basis  of  a most  simplified  fashion, 

e.g.  the  probability  of  error  for  the  binary  DPSK  system  has  always  been 

12  2 

based  on  the  expression  Pg  = ^ exp(-H  ),  where  H is  the  carrier-to-noise 

power  ratio  at  the  sampling  instant. 
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The  demodulator  configuration  of  a OPSK  system  requires  a phase- 
detector  (multiplier)  whose  inputs  are  "direct"  and  "one-bit  delayed"  versions 
of  the  signal.  In  reality,  the  inputs  of  the  phase-detector  have  different 
signal  power  levels  due  to  such  effects  as  phase  error,  delay  line 
attenuation,  and  inter-symbol  interference.  Furthermore,  the  noises 
between  the  two  multiplier  inputs  are  necessarily  correlated  due  to  the 
band-limited  nature  of  any  real  system.  When  these  practical  aspects  are 
considered,  the  error  behavior  of  the  DPSK  system  cannot  be  assessed  only 
on  the  basis  of  the  error  probability  expression  stated  above. 

The  lack  of  adequate  error  behavior  expressions  for  the  DPSK 
system  under  non-ideal  conditions  has  been  a consequence  of  the  non-linear 
nature  of  the  DPSK  demodulator.  Thus,  the  noise  behavior  at  the  decision 
time  is  not  amenable  to  a relatively  simple  description  such  as  coherent 
PSK  system.  As  a result,  the  analysis  of  DPSK  system  has  always  been 
based  on  a most  simplified  ideal  assumption.  In  particular,  the  noise 
correlation  between  the  two  phase-detector  inputs  has  been  totally  ignored 
in  the  idealized  assumptions. 

In  this  report  we  have  treated  the  case  of  the  DPSK  system  over  a 
hard-limiting  satellite  channel  with  power  imbalance  and  correlated  noise 
at  the  phase  detector  (multiplier).  A surprising  result  of  our  analysis 
is  that  the  error  performance  is  independent  of  the  noise  correlation  if 
the  a priori  symbol  probabilities  are  equal.  As  has  been  remarked 
elsewhere  [1],  this  is  a departure  from  previous  beliefs. 
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2.  ANALYSIS  MODEL 

The  binary  satellite  communication  system  under  consideration  is  shown 
in  Figure  1.  Our  main  objective  is  to  determine  the  error  performance 
of  the  system  under  the  assumptions  that  ( SNR) ^ f (SNR)2  at  the  phase 
detector  input  and  the  noises  n^t)  and  n^(t)  are  statistically  dependent. 
The  condition  ( SNR ) ^ f (SNR)2  arises  from  a difference  in  signal  powers 
due,  for  example,  to  intersymbol  interference  or  delay  circuit  phase 
error.  The  noise  correlation  reflects  the  band-limited  nature  of  a 
practical  system. 


The  original  transmitted  signals  are  defined  as  Sj(t)  and  S2 ( t ) : 


Sj(t)  = /2P ” cos(b>t  - 6j)  (1) 

S2(t)  = /2P^  cos [u>( t - T)  - 02] 

= /2F  cos (ut  - e2)  (2) 

where  the  bit  duration  is  T;  the  carrier  frequency  is  assumed  to  be 
selected  such  that  tuT  = 2nk,  k integer;  and  is  the  power  received  at 
the  satellite.  The  index  i = 1 is  associated  with  the  present  symbol  and 
i = 2 is  associated  with  the  previous  symbol.  Bandpass  Gaussian  noise 
nu(t)  is  also  present  on  the  uplink  with  ojj  = E|n^(t)|  . The  uplink  SNR 
is  defined  as  R^  - Py/o^. 


At  the  receiver  the  signal  is  corrupted  both  by  passage  through  the 
hard  limiter  and  by  additive  noise  on  the  downlink.  The  inputs  to  the 
phase  detector  (multiplier)  are: 

(3) 

(4) 


Uj(t)  = Sj(t)  + nj(t) 


u2(t)  = s2(t)  + n2(t) 
where  s_(t),  i = 1,  2,  are  the  signals  as  corrupted  by  passage  through  the 
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FIGURE  1 BLOCK  DIAGRAM  OF  DPSK  SYSTEM  OVER  A 
HARD-LIMITING  SATELLITE  CHANNEL 
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limiter  in  the  presence  of  uplink  noise  and  n.(t),  i = 1,  2,  are  additive 
downlink  noise.  We  may  write 

Sj(t)  = /2Py  cos(ut  + $j)  (5) 

s?(t)  = cos (ut  + (6) 

where  P^  and  P^  are  the  carrier  powers  at  the  phase  detector  and  the 
signal  phases  and  are  identically  distributed  random  variables  with 
the  conditional  density  function  [2],  [3] 


% <«IV 


(7) 


where  b 


Rk 

= Ek  nr  r(l + 1}  iFi(l;  k+1;  • 

is  the  signal  to  noise  ratio  at  the  input  to  the  limiter, 


Ek 


1,  k = 0 

2,  k > 0 


and  jFj(a;b;z)  is  the  confluent  hypergeometric  function  (Kummer's 
function)  of  parameters  a and  b and  argument  z.  The  bandpass  Gaussian 
downlink  noises  are  expressed  in  the  form 

hj(t)  = Xj ( t ) cos  ut  + Yj (t)  sin  ut  (8) 

n2(t)  = Xj(t-T)  cos  [u(t-T)]  + Yjlt-T)  sin[«(t-T)] 

= X?(t)  cos  ut  + Yp(t)  sin  ut  (9) 

and  the  noise  correlation  is  defined  by 

E{n1(t)n2(t))  = E{n1(t)n1(t-T))  = o*p(T) 


(10) 
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where  p(T)  is  the  normal  wed  correlation  with  T equal  to  the  symbol 
duration  and  the  noise  power  is  given  by 

a2  $ Eln^(t)}  = E{X*(t)}  = E{ Y^(t)) . i =1,2.  (11) 

The  decision  variable  y(t)  at  time  t is  the  output  of  the  zonally 
low-pass  filtered  version  of 
x(t)  = uj(t)  x u2(t) 

= [/2Pj  cos (ut  - + Xj(t)  cos  ait  + Yj(t)  sin  ait] 

x [/2P  cos ( ut  - $-)  + X,(t)  cos  ut  + Y,(t)  sin  ut]. 


where 


y(t)  = /PlP2  cos(*r*2)  + N(t)  (13) 

N(t)  = |x1(t)x2(t)  + |Yj(t)Y2(t)  + /fy7  Xj(t)  cos  *2 

+ A y2  Yj ( t)  sin  4>2  + /P^/2  X?(t)  cos  + JP^JZ  Y2(t)  sin  * 
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3.  STATISTICS  OF  THE  DECISION  VARIABLE 


The  probability  density  function  of  the  decision  variable  conditioned 
on  the  phase  difference  between  the  two  inputs  to  the  multiplier  is  a 
special  case  of  the  general  result  obtained  by  Miller  and  Lee  [4]: 


fy(y;pU)  = \ exp  [-(hZ  + hZ)] 


where 


5 S i [(’-?Kr  a (0?kr 

ik-J 

“'LfcKr- 


, y 2 0 


h3  = + h2  + Zhlh2  cos 

h4  = JTTTl  (hl  + h2  " Zhih2  cos  ** ■ 

2 

= SNR  at  input  number  1 of  the  multiplier, 

2 

= SNR  at  input  number  2 of  the  multiplier, 

$ = = the  difference  in  phase  between  the  narrow-band 

signals  at  the  multiplier  input, 

p = noise  correlation. 


■ E • 


The  probability  density  function  of  <i>  can  be  obtained  from 

<30 

%(PlVe2)  = Jf*1(6+a|el)f*2(ale2)da  (lf 

where  f.  and  f are  given  by  (7).  The  evaluation*  of  (15)  is  straight- 
s'] v 2 

forward  but  somewhat  lengthy  (see  Appendix  A),  with  the  result  being 


TTorlFn  alternate  approach  not  requiring  explicit  evaluation  of  (15),  see 
Appendix  D. 
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o. 


P-12 

\2n  1 


fielere2) = 


[21r-|e-(e1-e2)|] 

g b2^27,-|g^ere2)lj  cos[k|e-(ere2)|] 
- 2^sin[k|e-(e1-e2)|] 


§ S 


q=U  n k -q 

-q  sin[q|B-(e1-e2)|] 


k sin[k|B-(e1-02)|] 


1 6-j -02-2tt<  BS0-| -02+2tt 

1 0,  01-02+2tt<B<”.  (1G) 

An  important  aspect  of  the  p.d.f.  given  in  (16)  is  that  f ( 0 1 *©2)  depends 
only  on  the  difference  e^-e^  whlch  1S  the  information-bearing  parameter 
in  a DPSK  system.  We  can  thus  write  f (Bp  ,02)  as  This 

density  is  shown  in  Figure  2. 

The  probability  density  function  of  the  decision  variable  y 
conditioned  on  the  transmitted  symbol  (©  ^ ~ e2 ) with  downlink  noise 
correlation  as  a parameter  can  be  found  from 


co 

fy(y>p|ere2)  = / Ty(y,p 1 4>=&)p B 1 6-j -62)dB  . (17) 

The  result  of  (17)  is 
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FIGURE  2 PROBABILITY  DENSITY  FUNCTION  OF  PHASE  DIFFERENCE  BtTWEEN 
CONSECUTIVE  SYMBOLS  AT  OUTPUT  OF  HARD  LIMITER 
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fy(y;p|e,-B2)  - ’ exp  [-  -J-|]  ^ in  nA  (hl+hj)] 

d 


1 f1+Pl  fh24-h2 

4vraV'i+h2 


)]"  Vn 


where 


Vm  n = / exp[Y  cosb][1+Z  cose]m[l-Z  cosp]nf  ( e | e, -e2)de 

* -c=  9 (19) 


_ 2h1h2 
Z = —n *- 

hT+h2 


The  transmitted  symbol  is  "mark"  or  "space"  depending  on  the  phase 
difference  | -©2 1 =0  or  l e-j "e2 1 =1T » respectively.  Thus  the  following  conditional 
p.d.f.'s  are  obtained: 

fy(y;p|space)  = fy(y;P| I e1-e2|=o)  (20a) 


fy(y;p |niark)  = f (y;p|  |e1-e2|=7i) 


(20b) 


We  evaluate  (see  Appendix  B)  the  integral  in  (19)  for  both  space  (ej-e2=0) 
and  mark  (|ej-e2|=n)  by  using  the  binomial  theorem  arl  an  expansion  for 


powers  of  the  cosine,  with  the  results: 
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" 2 s 
Vm,n(spaCc)  “ 7 S E ^ ([S-U(5)-2w]/2} 


k=0  s=0  p=0 


s u(s)  b.c 


* (-D 


k p+ufs’ 


where 


* fI2p+u(s)-k^  + ]2p+u(s)+k^Y^ 
*(s)  2Fl^"m’"s;n"s+1  ;_1) 


s-uls) 
m 2 


V (mark)  = E ( [s-u(s)-2p]/2l 

m,n  1 k=0  s=0  p=0  ' 1 


2 

}kcp+u(s  ] 


^2p+u(s)-k^Y^  + ^p+ufsj+k^^ 


‘Is)  2F 1 (-n,“s »m_s+l *-l ) 


, s even 
, s odd; 


(21a) 


(21b) 


In(*)  is  the  modified  Bessel  function  of  order  n;  and  2Fj(. is  the 
Gauss  hypergeometric  function. 
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4 . CONDI T IONAL  E RROR  PROBABILITIES 

In  order  to  determine  the  error  probability,  we  must  know  the 
threshold  to  which  the  decision  variable  y(t)  is  compared.  In  a binary 
DPSK  system  the  transmitted  phase  difference  between  consecutive  symbols, 

el'e2’  either  0 or  111  • The  threshold  should  be  set  so  that  correct 
decisions  are  made  if  noise  is  absent.  In  the  absence  of  uplink  noise 
the  limiter  does  not  corrupt  signal  phase,  i.e.  4^  = and  = e^.  Then 
in  the  absence  of  downlink  noise  the  decision  variable  becomes 

y(*)|  no  noise  * cos<*r*2>  " cos  <Ve2>  <??> 

and  the  information  is  contained  in  the  sign  of  y(t).  Thus  the  threshold 
is  set  at  zero. 

From  (20a)  and  (20b)  one  can  obtain  the  conditional  probabil ities 
of  error  from  the  expressions 

0 

P(e;p|space)  = Prob{y<0| e^e^O)  = f -f  (y;p|e1-e2=0)  dy  (23a) 

— CD  ~ 

and 

co 

P(e;p|mark)  = Prob{y>0|  I©1-62I=ti)  = £ fy(y;p|  l®1~e2|=n)  dy.  (23b) 

To  compute  the  error  probabilities  we  need  to  define  the  downlink 

2 

signal  to  noise  ratio  Rd>  In  a DPSK  system  this  is  conveniently  taken  as 
the  SNR  at  the  direct  channel  input  to  the  multiplier.  In  our  notation, 

Rj  = h?-  (?4) 

Also,  since  we  are  considering  the  case  of  SNR  imbalance  at  the  multiplier 
inputs,  it  is  convenient  to  define  an  SNR  difference  measure  by 

? hp  (SNR)- 

* = ^ = (snrt;-  (25) 
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5.  ERROR  PROBABILITY  EXPRESSIONS 


Carrying  out  the  integrations  indicated  in  (23a)  and  (23b)  (see  Appendix 

2 2 2 

C)  and  using  the  definitions  in  (24)  and  (25)  to  replace  hj  by  R^  and  h‘ 
by  X Rjj  we  obtain  the  following  rather  formidable-appearing  expressions 
for  the  conditional  error  probabilities: 

s-  ' 

„ 2 “ n 

P(e;p,x|space)  = *(V)  exP \ 52  5T  52  52  52 

1-p  n=0  m=0  k=0  s=0  p=0 

r r r_i  ^(s) 

k p+u(s)v  ' /n+m\/n\ 

™ln!  rw^rp'  - l - )(s) 


I2p+o(s)-k^_p2  J + ^2p+u(s)+k^ _p2  ^ 
2Fi  (-n,l 2Fi (-m,-s  ;n-s+l  ;-l ) 
[lFl(!ikt,i-Ru)J  (26a) 


4 


13 


rsj  c 


]n  writing  (26a)  and  (26b),  we  have  made  use  of  the  relation  (2z) !=r(2z+l  ) = 

?7  1 9 

2 r(z+^)r(z+l)//ir  to  simplify  the  coefficients  arising  from  in  (21a) 

and  (21b).  For  alternate  ways  of  writing  (26a)  and  (26b),  see  Appendix  E. 


Making  use  of  the  relation  I .( -z)  = (-l )J I . ( z) , j an  integer,  and 
making  a change  of  notation  in  (26b)  by  replacing  m by  n and  vice  versa, 
we  observe  the  interesting  synmetry  property 

P(e;p,X | space) =P(e;-p,x| mark) . (27) 

The  total  unconditional  error  probability  is  the  weighted  sum  of 
(26a)  and  (26b)  given  by 

P(e;p,X)=PsP(e;p,k|space)+PmP(e;p,x|mark)  (20) 

where  P^  and  P^  are  the  a priori  probabilities  of  space  ("0")  and  mark  ("1"). 
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6 . COMPUTATIONAL  RESULTS 


We  have  computed  the  total  error  probability  for  the  case  of  equal 
a priori  symbol  probabilities,  i.e.  Pm=Ps=l/2.  In  doing  the  computations, 
the  direct  numerical  computation  of  (26a)  and  (26b)  proved  to  be  quite 
difficult  due  to  the  manner  in  which  the  multiple  series  converge  and  the 
effects  of  finite  word  length  in  the  computer.  We  found  that  numerical 
integration  of  (19)  in  conjunction  with  (18),  (23a)  and  (23b)  proved  to 
be  a more  computationally  efficient  scheme  for  obtaining  numerical  results, 
and  this  was  the  procedure  used  in  our  computations.  The  program  listing 
is  contained  in  Appendix  G. 

Our  numerical  evaluations  of  the  total  error  probability  were  performed 

2 

for  several  values  of  power  imbalance*  X and  noise  correlation  p.  We  have 

observed  the  result  that  the  computed  total  error  probabilities  are  inde- 

2 

pendent  of  p for  all  values  of  X considered.  However,  the  conditional 
error  probabilities  were  dependent  upon  the  noise  correlation.  When  these 
conditional  error  probabilities  were  added  with  equal  weights  (Pm=Ps=l/2), 
the  results  coincided  with  the  value  of  the  error  probability  for  p=0,  as 
shown  in  Figure  3.  A similar  numerical  result  was  obtained  for  a terrestrial 
link  [5]  (see  also  Appendices  H and  I).  As  in  the  case  of  the  terrestrial 
link,  the  mathematical  complexity  prevented  analytical  verification  of  this 
result.  It  should  be  stated  here  that  when  one  observes  the  conditional 
probabilities  of  error  as  given  by  (26a)  and  (26b),  it  appears  certain  that 
the  total  unconditional  probability  of  error  should  depend  on  noise  correlation 
p.  In  fact,  it  has  been  remarked  in  the  previous  publications  [6],  [7]  that 
the  error  probability  of  a binary  DPSK  system  depends  upon  the  correlation 


*The  power  imbalance  can  arise  from  intersymbol  interference  and  other  effects, 
as  discussed  in  Appendix  F.  15 
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UPLINK  SNR  R*  (dB) 

FIGURE  3 INFLUENCE  OF  NOISE  CORRELATION  ON  CONDITIONAL  ERROR  PROBABILITIES 
FOR  POWER  IMBALANCE  X2  = -1.5  dB  AND  DOWNLINK  SNR  R?  = 7 dB 
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if  there  is  intersymbol  interference  (a  case  of  SNR  imbalance).  The  remarks, 
it  should  be  pointed  out,  were  based  on  the  observations  of  the  equations 
rather  than  on  the  proofs. 

Figure  4 shows  the  total  error  probability  versus  uplink  SNR  R^  for 

several  values  of  \2  and  R^.  For  comparison,  the  case  of  infinite  downlink 

SNR  (Rd=»)  is  also  plotted;  this  is  identical  to  the  terrestrial  link  with 

2 

no  power  imbalance  (x  =0  dB). 

The  curves  shown  in  Figure  4 corresponding  to  x2=0  dB  for  each  case 

2 

of  Rd  (downlink  SNR)  represent  the  ideal  case  of  no  power  imbalance 
at  the  phase  detector  input.  Our  curves  for  these  special  (ideal)  cases 
are  identical  to  the  results  reported  by  Weinberg  [8]. 
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7.  CONCLUSIONS 

In  this  report  we  have  presented  the  error  behavior  of  a binary  DPSK 
system  over  a hard-limiting  satellite  channel  under  the  influence  of  an 
SNR  imbalance  at  the  phase  detector  (multiplier)  inputs  and  noise  correlation. 
The  graphically  presented  curves  for  error  probability  are  applicable  to 
evaluation  of  system  performance  when  the  SNR  imbalance  at  the  phase  detector 
is  known. 

Our  numerical  results  show  that  the  performance  of  binary  DPSK  over 
the  hard-limiting  satellite  channel  does  not  depend  upon  noise  correlation 
when  the  a prion  symbol  probabilities  are  equal,  regardless  of  SNR 
imbalance.  The  noise  correlation  has  an  effect  only  when  the  symbol 
probabilities  are  unequal. 

The  mathematical  complexity  of  the  error  rate  expressions  has,  thus 
far,  prevented  an  analytical  verification  of  the  numerical  results  that 
the  noise  correlation  has  no  effect  on  performance  when  the  symbol 
probabilities  are  equal.  The  mathematical  formulation  for  the  error  rate, 
however,  is  possibly  amenable  to  further  investigation.  Appendix  E points 
out  a few  first  steps  which  may  lead  to  further  analytical  investigation 
of  the  error  rate  properties;  however,  the  mathematical  relations  involved 
are  in  an  area  of  mathematics,  e.g.  generalized  multiple  hypergeometric 
functions,  which  is  not  yet  fully  developed.  Further  investigations  may 
be  mathematically  interesting,  but  the  directions  of  such  investigations 
and  the  immediacy  of  practical  results  of  such  investigations  are  by  no 
means  clearly  evident  at  this  time. 


J.  S.  LEE  ASSOCIATES,  INC 


APPENDIX  A 

DERIVATION  OF  PROBABILITY  DENSITY  FUNCTION  OF  * = 4^  - 


The  pdf  of  i=l,2,  is  given  by  equation  (7)  of  the  main  text. 
We  must  evaluate 


f^(e|ere2)  -r  ^ (e+alei  >f (a I e2 >dc 


(A-l) 


In  view  of  the  restricted  range  over  which  the  f 's  are  non-zero,  (A-l) 

$i 

must  be  treated  as  four  separate  cases  (see  Figure  A-l).  If  -3+e^+7i<e2-7i, 

or  0.-e9+2ir<e,  then  there  is  no  overlap  of  the  two  f 's  and  f =0.  In  the 
\ e.  9i  4 

second  case,  we  have  e^-iK-e+e^+TKe^+n,  which  implies  that 

/•e.+n-p 

f*(B|e.,e2)  = / 1 f.  (e+ale^f.  (a|82)da,  e1-e2<B<e1-e2+21,. 

V*  (A-2) 

In  the  third  case,  -p+Oj -u<^0?+n< -B+Oj+ti  , which  implies  that 

f.Uh.e,,)  = f02  ” % (B+ale^f.  (a|e2)da,  ej-e^nsBsej-eg  . 

el-e  (A-3) 

In  the  fourth  case  -B+8^-7i>e2+7i , or  B<0j-e2-2ji,  there  is  no  overlap  and 

f ♦ ( b I ei  ,e2)  = 0. 

Substituting  (7)  into  (A-3)  and  (A-2),  we  find 


UO  CO 


f,(B|e..e,)  ■ (i-j  f2*’  cos[k(B<»-e,)]cos[q(a-e2)]d., 

9 Id  \di\  / k=o  q=Q  K q -,6j-ii-B 


6j-e2-2-n<Bs0j-e2 

(A-4) 


A-l 
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B+0J- 


P( a ) p(a) 


c)  Case  3 


d)  Case  4 


FIGURE  A-l  DETERMINATION  OF  INTEGRATION  LIMITS 
FOR  PROBABILITY  DENSITY  OF  $ AS 
A FUNCTION  OF  B 1 1 
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2 ^ W y-0-j  +*n-  P 

%(ple]  »e2)  = (2^)  EE  bkbq  J cos[k(8+a-01  )]cos[q(a-e2)]da, 


k=0  q=0  e2-if 


e1-e2je^e1-e2+2Ti  (A-5) 


where 


? r(t+  ])  1 Fi (2* k+1;  _Ru)- 


bk  = Ek  k 


Now  consider  the  integral  in  (A-4): 


up ' n 

I = J COS  [k ( 0+a-0 -j  )J  cos[q(a-02)]  da. 


(A-6) 


e-j-n-B 


We  have  three  cases,  determined  by  the  parameters  k and  q.  When  k=q=0, 
the  integrand  of  (A-6)  becomes  equal  to  one,  and  we  have 


(A-  7 ) 


I = 2n+B+02-e^  > k=q=0. 

The  second  case  is  k=q^0,  for  which  we  apply  a trigonometric 
identity  to  write  (A-6)  as 


e2+Tt  e2+u 

I = cos  ( ke+ke2-ke1 ) da  + j / cos(2ka+kB-ke1-ke2)  da. 


e -j  - 7t-b 


Integrating  (A-8)  we  obtain 


6^ -t-6 


(A-8) 


A-3 
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1 = 


2 cos  (kH+kO^-ko^ ) «n  + p) 

+ L [sin(2ke2+2kn+k3-k0^-ko2)  - sin(2kG^-2ku-2kp*kp-kOj-kez)] 
2 cos  [k^k(o?-o1)]  (2n+p+o?-e1 ) 


[sin(kp+ko^-ke^ ) + sin(kp+ko2-kPj )] 


(A-9a) 


which  can  be  written  as 
1 


1 = ^(2ii46+e2-e1)cos[ke+k(o2-e1)]  + sSnCke+k^o^-e, ) 3, 

k=qyo.  (A-9b) 

In  the  third  case,  we  have  k/q.  Using  a trigonometric  identity, 

(A-6)  can  be  written  as 

i r 2*'  1 fc" 

I ~ 2 I cos[(k-q)a+kp+q02-k0j ]d«  + j / cost  (k+q )u+kp-kP j -qe,,]da 


Integrating  (A- 10)  yields 
1 


6j -n-B 


(A- 10) 


I = 2“(k-q^ {S1*  n[ke2-q02+(k-q)Ti+ke+qO2-ke1  ] 

- sin[k0.-q0^  - (k-q^-kP+qP+kB+qP^kO-i  ] ) 


+ ?(k]+qy  lsin[k©24qe2+(k+q)n4kB-ke1-qe2] 

- sin[k6^+qP^- (k+qJn-kP-qB+kP-kO^-qP^]} 

= 2(k-q)  <sin[(k-q)it+ ke +k(P2-ei)]  - sin[-(k-q)i'+qP+q(02-01)]} 

+ ?(k+qy  {SH\[(k4q)n4kP+k(P2-ei)3  - sint-lk+q^-qP-q^-Pj )]) 


A-4 


which  simplifies  to 
1 


I = 2(k-qT  {sint  (k-q)T.  + kP+k(02-Oi  )3  + sin[ ( k -q ) tt -q p-q ( e2~ e 1 )]} 

+ 2(kW  {sint(k+cl)Tr+ke+k(02'ei  >3  + sin[(k+q)1>4qe+q(e2-e1 )]}. 

(A-ll) 

Making  use  of  the  identity  sin(x+y)  = sin  x cos  y + cos  x sin  y,  and 

1/ 

recognizing  that  cos  (K")  = (-1)  , (A-ll)  can  be  written  as 

, = i-Uq — sin[(k-^)e  + (k^ ) ( e2"  ei )]  cos  [(k|q)e  + (k2S)(e2-0i )] 

+ 4+q “ Sin  [(^)6+  (4a)(V0l)]  C0S  [(k2q)B  + (lfl)(V6l)]- 

(A- 12) 

Noting  that  (-l)k~q  = (-l)kq+^q  = (-l)k+q  and  that  2sinxcosy  = 
sin  + sin^-^j,  (A-12)  can  be  written  in  the  form 
f-llk+q 

I = A~g ' £ (k  sin[ke+k(o2-e1 )]  -q  sin[q8+q(e2-e1 )]} , k/q.  (A- 13) 

We  now  substitute  ( A- 7 ) » ( A- 9 ) , and  (A- 1 3)  into  ( A-4 ) to  obtain, 
replacing  (e2~ 0-, ) by  -(e^o^). 


Vb|V82>  [2"+B-(e1-e2)] 


(if  E 


r^TT+e-  (e-j  - e2 ) 


k=l 


] 


cosOe-k^-e^)] 


sinCk6-k(e^-e2)3^  + (2V)  z±t\\y 

' u = n n-r\  * “ 


k+q 

r 


k=0  q=0 

k/q 


• {k  sin  [ke-k(e1 -e^)]  - q sin  [q8-q(e1-e2)]l, 


8j-e2-2n<B<ej-e2. 


(A- 14) 


A-5 
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In  ( A- 1 4 ) the  first  term  in  square  brackets  accounts  for  the  k=q=0  term 
of  (A-4),  the  single  summation  accounts  for  the  main  diagonal  k=q,  and 
the  double  sum  encompasses  the  remaining  terms  of  (A-4). 


We  observe  from  (A- 13)  that  for  q^k,  I ( k=kQ  , q-qQ ) = I (k=qQ,  q=kQ) • 


Thus  the  double  sum  in  (A-14)  can  be  reduced  to  one  infinite  sum  and 
one  finite  sum: 


V-'  M)k+0 


bkbq  ,2  ~2~  {k  sin  tke'k(0re2^  -q  sin  CqB-q(e1-e2)]}. 


k=l  q=0 
Thus  (A-14)  becomes 


kt_  & 

-q 


Ve|el,e2)  = (i)  [2«*B-(ere2)] 


(lft<  jp^.] 


cos[kB-k(e^-e2)] 


k=l 


+ i sin[ke-k(e1-e2)] 


k-1 


+ 2(s)  S S bkbq  ,k  s1n  l*e-k(ere2)] 

k=l  n=Q  k "q 


- q sin  [qe-q (e] -e2)3 > , 


e^-e2-2-n<B$e^-e2. 


(A- 15) 


From  (A-5),  we  need  to  evaluate 

J = J cos(ko+kB-ke^ ) cos(qa-qe,,)  da 


(A- 16) 


02_ir 


which  differs  from  (A-6)  only  in  the  limits  of  integration.  Again 
there  are  three  cases  to  consider. 


A-6 


k 


. .• 

I 


l: 
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If  k=q=0,  then  ( A- 1 6 ) becomes 


J = 


■/ 


+n-0 


da,  k=q=0 


( A- 17) 


&2~v 


which  yields 


J = 2n-B-(e2  - 6} )»  k=q=0. 


( A- 18) 


If  k=q^0,  then  (A- 16)  becomes 


J = 


/ 


,+ti-B 


cos(ka+k0-k6j ) COS (ka-ke^)  da. 


(A- 19) 


02“TI 


Using  the  identity  cos x cos y = ^ [cos (x+y ) + cos (x-y)]  we  can  write 
(A- 19)  as 


= */ 


0^+7T-0 


cos[kB+k(0^-02)]  da  + 


02~it 


u 


0,+n-  B 


cos[2ka+kB-k02-ke1  ] da. 


02_n 


(A- 20 ) 


Integrating  (A-20), 

J = ^ [2r-B~ ( 6 2_ e -j  ) 3 cosCkB+kCe^e^] 

+ {sin[2ke1+2kTi-2kB+kB-k02-ke1] 
- sin[2ke2-2kTi+kB-ke2-ke1]) 
2tt-b-  (e2-e1 ) 


cosCks+k^-e-j)] 


+ ^L  (si n[-kB-k(e2-e1 ) 3 - sin[kB+k(e2*e1 )]} 


(A-21 ) 


1'! 


■i 


1 ■ 


> i 


. * J: 
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which  simplifies  to 


J = 1 [2„-e+(ere2)]  cos[ke-k(ere2)]  - -L  sintkB-k^-e^],  k=q^0. 


( A- 22 ) 


If  Mq,  then  (A- 16)  becomes 


0 j+n-8  u.-j  ■ » _ p 

= \j  cos[(k-q)a+kB+q02‘k6^  ] da  + ^ J cos[(k+q)a+k0-k0^-q02]  da 


0,  +TI-0 


(A-23) 


which  integrates  to 


J = 2fk-qT  ^s’ntkei_<^0i  + tk_t^)7I_ke+^f5+kp+c^e2-,cel^ 

- sin[k02-q02-(k-q)Ti  + kB+q02-k0^]} 

+ 2"(kVq7  {sin[k01+qe1  + (k+q)n-kB-qB+kB-k01-q02] 

- sin[k@2+q02- (k+q)r+kB-k0^-q02]l 

= {sin[(k-q)rr+qB+q(e2-01 )]  - sinf-fk-qMkB+k^-Oj )]} 

+ 2(kVqT  ^sin[(k+cl)11't1B-cl(e2'el  ^ ‘ sin[-(k+q)Ti+kB+k(e2-01 )]} 

= 2Xk‘-qT  {sint (k-qjT’+qe+q (e2_0l  ^ + sin[(k-q)Ti-k6-k(e2-01  )D> 

+ -2(k+q-y  { si n[ (k+q )Ti-qB-q (e2-e1 )]  + s i n [ ( k+q ) tt - k b~ k ( 02~ 0 1 )]} 

= 2(k~-qT  {cos^k"q^^  sin[qB+q(e2-e1)]  - cos[(k-q)*]  s i n [kB+ k ( ©2~ 0 1 )]} 

+ 2(k+qJ  {-cos[(k+q),]  sintqB+qfe^Oj )]  - cos[(k+q)n]  s i n [kp+k ( e2~ e ] )U) 

(A-24) 


to  introduce  a finite  upper  limit  on  the  inner  summation. 

; f 

In  (A- 15),  we  can  rewrite  the  limits  of  applicability  as 

-2TT<e-(e1-e2)<o 

' l 

and  in  (A- 26 ) we  can  rewrite  the  limits  of  applicability  as 

0jB-(e1-e2)$2n. 

Thus  ( A- 1 5 ) and  (A- 26 ) apply  to  adjacent  regions,  depending  upon 
whether  B-(0^-02)<O  or  B-(0j-6p)>O,  respectively.  Letting  if>  = e- (©^ -e^) 
in  (A- 15)  and  (A-26),  we  can  see  that  ( A- 1 5 ) is  merely  (A- 26 ) with 
♦ replaced  by  -<K  Since  (A- 15)  applies  for  over  the  range 

of  applicability.  For  ^>0,  |i|<|=i|<  and  thus  ( A- 1 5 ) and  (A-26)  can  be 
combined  to  give 


A-9 
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o. 


£)  j[2.-|B-(ere2)|] 


f(e | e1 »e2)  = 


Eoir^n- | e- ( o, - ) n 

bk|L 2 J costklp-(°re?) n 

k=l 

- 2k  sin  [k | g- (©i-e^) | ] | 


k-1 


-l  Yh H bkbqtT  2- jksin[k|e-(ere2)|] 


k=  1 q=0 


\o. 


k -q^  / 

-q  sin  [q|e-(e1-e2)|]||, 

0i  - 02*2ti<b$  6i  ■( 

0-|-0j+2nf  g<". 


2+2t 


A-10 
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APPENDIX  B 

EVALUATION  OF  EQUATION  (19) 

We  will  evaluate  (19)  for  two  cases:  I ©j -©Pl  s 0 ("space")  and 

[Oj-e^  = v ("mark").  First  we  treat  the  case  of  "space"  being  trans- 
mitted. 

B.i  "space"  (|ere2|  = 0) 

We  use  the  binomial  theorem  twice  to  write  (19)  of  the  main  text 


CT1 

)>n  * J exp[Ycose]  f ^ ( e 1 0 ) 


[E  oiwj£  c 


g(-Z)rcosrejdB 
r=0  J(B-1) 


where  we  have  taken  into  account  the  finite  range  over  which  f^eJO) 

does  not  vanish.  Interchanging  the  order  of  summation  and  integration, 

which  is  clearly  permissable  due  to  the  finite  limits, 
m n 2v 


f=0  r=0 


(l)  (r)( -1 ) P 2i+rJ exp [Ycosb]  cos£+r6  0>d&'.  (B-2) 


Denoting  the  integral  in  (B-2)  as  Ws(£+r),  using  (16)  of  the  main  text 

for  f (e|0),  and  splitting  the  range  of  integration  into  two  subintervals 
v 

for  positive  and  negative  B to  enable  us  to  eliminate  the  absolute 
value  signs  from  the  integrands,  we. have 


WS(M)  = f cosMBeYcosB^r)2|(21i+B)- 

-2ti  ' 

+ cos(ke)  4^sin(kB)l 

i.  _ n J 


B-l 
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k=l  q=0 


b.b  (-1 
k q — ? 


(-l)k+q  T l) 

— ? k sin(kB)  - q sin(qe)  > ds 
k - q L J) 


Pv 

♦ / ’ cosV“sB(^2{(2.-6) 

, b2  [(  2*jj8  Ws(!<e)  - sin(kB) 
k=l  *■ 


k=l  q=0 


bkbq  — ? [k  sin(ke)  " q sin(qB)j|d6  (B-3) 


where  we  have  defined  M=£+r  for  notational  simplicity.  If  we  make 
the  change  of  variable  y = -B  in  the  first  integral  of  (B-3)  we  find 
that  it  equals  the  second  integral.  Interchanging  the  order  of  summation 
and  integration,  we  have: 


WS(M>*  2"COSM3eYcosBdB  'f 


B cosMBeYcOSBdB 


°1^  r-  /*2ti 

^ b^  J o cos  kecosMe  eYcosede  - J b COS  kBcosMgeYcosBdB 

if  'sink6cosMBeVcos6d6j 

2Sj3bkbq^7’[kl  sinkBCOSH8eVcOS8dB 

k=l  0=0  K -q  L 0 


q / sinqBcosM6eYcosB 
J0 


de]|- 


(B-4) 
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In  (B-4)  we  make  the  change  of  variable  a=B-n  and  make  use  of  the  identity 


COs(Qa  + Qn)  = (-l)^COSQa 
sin(Qu  + Qn)  = (-l)QsinQa 


to  write 


Q integer 


(B-5) 


VM>-  Hhf  {*•<-»"£  nMcosVVcosv 

-n  f (-l)MCOSMae'YC°Sada 

J-  tl 


+ [ * f (-l)kCOSka(-l)McosMae'YCOSada 

W 1 ^ 

a(-l)kCOSka(-l)MCOSMae"YCOSada 

r 


1 f , , ,k  , , , M -Ycosa. 
-if  n( - 1 ) coska(-l)  COS  ae  da 

n 

/.  tt 

(-1)  sinka(-l)  COS  ae  da 

n 


*'[>£ 


k=l  q=0 


/ M -Ycosa. 

(-1)  sinka(-l)  COS  ae  da 


q f (-l)qsinqa(-l)McosMae  Yc0SndaJ|. 


(B-6) 


In  (B-6)  the  second,  fifth,  seventh,  eighth,  and  ninth  integrals  vanish  because 
the  integrands  are  odd  functions  and  the  limits  are  symmetric  about  zero.  Also, 
we  may  write  the  integrals  of  even  functions  as  twice  the  integrals  over 
(0,n)  to  give 
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;<«>  ■ «(£)*  «-»)"  I- / 


iMuda  (-l)k 


n / -Ycosa  M . , 

• j I e cos  aCOSkada\. 


(B-7) 


Combining  formulas  for  the  even  and  odd  powers  of  the  cosine  [9,1.320.5 
and  1.320.7]  we  can  expand  the  power  of  the  cosine  in  terms  of  cosines 


of  multiple  arguments: 


M-u(M) 

r 


m i y* 

COS  c = W Z-J 
2 p=0 


'u-Hj(M) 


([M-u(M)-2u]/2^ 


cos{  L 2p  ( M ) D c } 


(B-8) 


. ■{*•»■ 
u {2, 

(0,  H 

u(H)  = |ls  M 


odd  . 


Using  the  identity 


cos  z cos  c = | [cos(zn)  + cos(z-c)] 
in  connection  with  (B-8)  allows  us  to  write  (B-7)  as 


(B-9) 


M-u(M) 

2 


W$(M)  = 4 


GJ  (_1)M  * f 6 YC°Sa  J M Cu  ([M-u  00-21.1/?)  C0S{[2w+u(M)]a}da 

JO  ‘ u=0  V ' 


M-u(M) 

2 


+ 2bk  (_1)  2 [ e S 2M  cp+u(M)  f[M-u(M)-2u]/2^ 

k=l  Jo  * t^O  V ' 
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•[|  cos{[2u+u(M)+k]a)  + ^ cos{  [2p+v(M)-k]a}j 

Interchanging  the  order  of  summation  and  integration  in  (B-10),  we 
obtain 


(B-10) 


M-o(M) 

2 


W$(M)  = A 


(?1)  X)  Cm  + u(M)  ^[M-u(M)-2y  ]/^ 


mj  e acos{ [2y-hj(M) ]a}da 
JO 


M-u(M) 


— ■ i 

*X/  bk  2 pi  X/  EP  +u(M)  Z[M-u(M)-2u]/2) 

k=l  y = 0 V / 


7! 

ij 

7 

**/ 

Jo 


e Y °Sacos{ [2p+u(M)+k]a}da 


e'YC0Sac0S{[2p-K)(M)-k]a}da 


( B - 1 1 ) 


All  of  the  integrals  remaining  in  (B-ll)  are  of  the  form  [10,  9.16.9] 

(B-12) 


l 


g-Ycoso  cosNa  da  = ^ (_Y). 

N 


Also  noting  that  IN(rY)  = ( - 1 ) N I N ( Y ) for  N integer,  (B-ll)  becomes 


V"i  ■ < (s)  <-1>M 


M-u(M) 

O 2 

JL.  / nu(M) 

„M 


X!  \ +u(M)  ^’M-u(r1)-2y 3/2)  12jj+o (M) 

u=0  ' / 


(Y 


M-u(M) 

vA 


+ XI  bk('1)k  2 2M+l  L % +u(H)  ([M-u(M)-2w]/2) 

k=l  y=0  ' / 

• l(-i)2p+u^+ki  , , (y)  + (v\l 

V 11  ^y+vM+k1^  1 I2p-K*(H)-kCYy 


B-5 


(B-13) 
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Since  ( -1 ) ^ = (-1)  ^ and  (-l)2w  = 1,  ( B - 1 3 ) becomes 


M-u(M) 


W$(M) 


= A 


(if  ( 1)n(»?(-i)“(")'sS  t / H \ 

\2ti/  v ; cp +o(M)  ^[M-u(M)-2p]/2/ 

U=0 


2v+v (M) 


M-u(M) 


+ ^ b2  / M v 

L k ?M+2  Z-,  ep  + u(M)  (LM-u(M)-2m]/2 

k = l C g=0  ' 

[I2^+v(M)+k(Y)  + ^uMMj-k^  | • 


( B - 1 4 ) 


Noting  that  b^=  1,  the  single  summation  in  (B-1A)  can  be  brought  inside 
the  double  summation  by  introducing  the  factor  1/c^  since  for  k=0  the 
two  modified  Bessel  functions  are  identical.  Thus, 


WS(M) 


(-1) 


M+u(M) 

7M+1 


M-u(M) 
2 


E E El.  %+u(M)  bk  ('M-u(M)-2u]/2) 

k=0  u=0  K ' / 


[I2u+u(M)-k(V)  + I2y+u(M)+ktY,j  • 

Putting  ( B -1 5 ) into  (B-2)  we  arrive  at  the  form 
■C-Hr-uU+r) 


( B - 1 5 ) 


m n °° 


m 


.»=EE  E E G) 

1=0  r=0  k=0  v= 0 


h2 

. k cu+u(l+r) 
,f+r+l 
2 ek 


[ l'r  ) 
\C£+r-u(£+r)-2u]/2/ 

* [!2y+uU+r)-k(Y^  + 1 2P ( £+r ) +k { Y * 


( B - 1 6 ) 


(Y) 
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Noting  that 


(b)  = 0 1 f b< 0 or  b> a , 


( B - 17 ) 


we  can  extend  the  limits  of  the  summations  over  £.,  r,  and  pin  (B-16) 
to  infinity.  Then  we  can  interchange  the  order  of  summations  and  use 
Bailey's  theorem  [11, pp.  58-59]  to  replace  the  doubly  infinite  sums 
over  l and  r with  one  infinite  sum  and  one  finite  sum.  Thus,  with 
s = f+r  and  t=£,  (B-16)  becomes 


i,n  = \ EE  EE  5Z(rs-p(s)-2p]/2)  (2 
k=0  p=0  s=0 


A (_i)u(s)  b2  Cvi 


k e, 


*2p  +u(s)-k^ 


+ I2„*u(s)+k<y>l  • <B-18) 

J t=0. 


It  can  be  shown  [12, p.  17]  that  the  finite  sum  over  t in  (B-18)  is  a 
Gauss  hypergeometric  function  with  a negative  numerator  parameter. 


E ft)  (s-t)*-"1  ■ (s) 


iB-19) 


Putting  (B-19)  into  (B-18)  and  again  using  (B-17)  to  introduce  finite 
limits  in  the  transformed  summations,  we  have  finally  that 

s-u(s) 

°°  n 2 2 

Vm  n(space)  = j EE  E ([s-Ws)5-?^)  (f)S(-Du(s)  ^-£y+u(s' 

k=0  s=0  p=0  ' ’ Ek 


I 2p>u(s) 


-k(Y)  + I2p+u(s)+k(Y)](s)2Fl{',n*'s;n"s+1;'1)  (3‘ 


which  is  (21a)  of  the  main  text. 
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B.2  “MARK"  ( I e e2|  = Tf) 

Again  we  use  the  binomial  theorem  twice  and  interchange  the  order 

of  integration  and  summation  to  obtain,  for  "mark", 

m n e1-e2+2n 

' fni\/n\,  ^r.£4r  C T 1 . 

cos  8 f^V8||e1-e2J=7r)dB. 

(B-21) 


m n 12 

Vm,n  = 5^  X/  (I'Xr)^'1^2  ^ f exp[ycose]  — '&+r- 

1=0  r=0  J 


0 j-02~27r 


Equation  (B-21)  is  identical  to  (B-2)  except  that  |8j-e2|=  n . 

Denote  the  integral  in  (B-21)  as  W^+r).  Since  |ere2|-  , „e  have 

two  cases,  namely  Oj-Oge  -TI  an(*  ei'e2  = n * which  occur  with  equal 
probability.  Therefore, 


7T 

WM^r)  = \ / e^p[YcosB|cos4+rB  f(B|e, 

^ J * 1 


^(b | e2-e2=-  •rr)dB 


1 f 

+ \ J exp  Jycosej  cos  ^+rB  f ^,(  B J e j - 

-TT 


92=n)dB. 


(B -22) 


Using  (16)  of  the  main  text  for  ^(bJOj-  e2),  letting  M = £+r, 
can  write  for  the  first  integral  in  (B-2Z) 


we 


(-^-)  f eYcOSBcosMB  !(3ti+b)  + ^2  b^  j-^6-  cos[k(B+*3 

J- 3tt  ' k^l  1 


+ i sin[k(B+Ti)} 


k-1 


k=l  q=-0  q k - q 


|k  sin[k(e+Trj]  - q sin  [q(e+t,j]  jdB  + (^-) 


eYcosBcosM  B 


(B-23) 
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{CO 

(*-®)  + bk  | (^T-)  cosjk(B+n)j  - sin  jk(e+«)j| 

CD  (("I 

E E v,  ^ 7 jk  sin  jk(B+*)j  - P sinjq(p+T!)j|| 


-2 


and  for  the  second  integral  in  (B-22) 


de 


f eYcOS6cosM6 |(n+6,  + £^j(^)cos[k(6.,)]^sin[k(6.^j 

-»  1 k=l 

® k-1 

-EE 


-71 

® k-1 

k«l  q=0  bkbc>  s14(8-’>]  - <5  sin[q(6-,)]|jde 

r3n 


2 / o'  . 

+ (^)  J eYcosBcos^|(3TT.6)  + ^ bk  cos[k^-i>)] 

k **1 

v °°  k-1 

- i sfn[k(8-,)]|-2 

1 k=l  q=0 


k 'I  I? - q2 


k sinJk(B-n)) - q sin  jq(fi-n  jjjjdB. 


(B-24) 


If  the  change  of  variable  y=-B  is  made  in  (B-23)  we  find  that 
W„(M)  = W*(M).  Thus  the  equally  weighted  sum  of  the  two  integrals  in 
(B-22)  equals  the  unweighted  value  of  either  integral. 

Using  the  form  in  (B-24)  for  Wm(M),  we  make  the  change  of 
variable  y=B-h.  Since  cos(y+h)  = - cosy,  we  have 


B-9 
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o 

VM )=\^r)/  e Yco^(-iAosmy|(2ti+y}  +£  bk[(“112JL)C0S^kY^  sin(k^^] 

\ 'J-2*  ' k=l 

+ 2 E 2 b,bj4i *-  |"k  sin(ky)  - q sin(qy)|/ dy 

k=l  q=0  K(>r-q  L ->) 

♦ (h)7  eVC°SY(-lAoA  |(2-y)  ♦ E b^(^)cos(k,)-^  Sin(kr)] 


°°  k-^  k+q  r ) 

- 2 2 S bkbq^2  _ 2 [K  sin<kY>  - Q sin(qy)  J dY.  (B-25) 

Comparing  (B-25)  with  (B-3)  we  see  that  the  only  differences  are  a 

M 

multiplicative  factor  of  (-1)  and  a change  of  sign  in  the  argument 
of  the  exponential.  This  latter  difference  will  lead  to  a change  in 
sign  of  the  arguments  of  the  modified  Bessel  functions,  allowing  us  to 
write  the  results  of  (B-25)  by  analogy  to  (B-15)  as 


, ,,2M+u(M) 

V">= 


M-u(M) 

OP  ? 

Z E 

k=0  u=0 


'wtv(M)  ^2 
b. 


([M-u(H)-2u  ]/zJ 


{I2P+u(M)-k("Y)  + ^p+offO+k^j 


(B-26) 


Since  ln(-x)  = (-1)  ln(x).  for  n an  integer. 


W„(H)  = M 


'21“  D £<^^2/  \ 

n 1 k=0  p=0  ck  k\  / 


M-u(M) 

2 


[<-1)'k,2ukv(M).k1',>  1 >-»k'2l,*v(«)kk<''>]  . 


(B-27) 
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Furthermore,  since  an  even  power  of  -1  is  equal  to  1,  and  ( -1 ) ' 


(-1)  , (B-27)  becomes 


M-u(H) 
2 


wm(m) 


,M+ 1 


V V VLylM)  (.i)k  2/  M \ 

U U \ bk(rM-u(M)-2p]/2j 


+ I?u+u(M)+k^Y^]  • 


( B -28) 


Therefore,  from  (B-21)  and  (B-28)  when  a "mark"  is  transmitted 


m n 


m,n 


- E E (X)<-»r  z 


l+r 


1+r-vU+r) 

L . V V Zv+»U+r) 
,£+r+l  JLt  e," 


1=0  r=0 


k=0  v=0 


(-nk 


/ l+r 

([£+r-u(.£+r)-2p]/2 


?) 


[I2v+u(/+r)-k^V^  + ^p+uU+rJ+k^H  * 


(B-29) 


We  again  use  (B-17)  to  introduce  infinite  limits  on  the  summations  and 
interchange  the  order  of  summations,  summing  first  over  k and  p then 
over  l and  r.  Then  we  apply  Bailey's  theorem  to  diagonalize  the  sums 
over  l and  r.  Letting  s =£+r  and  t=r,  (B-29)  becomes 

Vn ■ IE  E E E 

k=0  p=0  s=0  t=0 

*([s-u(s)-2p]/^  [I2p+u(s)-k(Y)  + I2vl+u(s)+k^Y|]  * 


(B-30) 


1 


The  summation  over  t in  (B-30)  becomes  a Gauss  hypergeometric  function 
[I2,p,17]  as  was  shown  in  (B-19).  Therefore, 


_ m s 

•m,n  " 2 E E E L(,UJ  (f)  (-Dk  b l 

k=0  s=0  p=0  ' ’ k 


[I2u+u{s)-k^V^  + I2p+u(s)+k^Yi| 


(s)  2Fl^"n,‘s;m's+1;'1^ 


(B-31) 


where  we  have  interchanged  the  order  of  summations  over  s and  v and 
have  again  used  (B-17)  to  introduce  finite  summation  limits. 


B - 1 2 
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APPENDIX  C 


EVALUATION  OF  EQUATIONS  (23a)  AND  (23b) 


C.l  EVALUATION  OF  (23a) 

Using  (20a)  of  the  main  text  in  (23a)  of  the  main  text, 

fe)  K-D] 

-Co  rl  »-  ■ M J a a 


n=0  m=0 


where  Vm>n  is  not  a function  of  y (see  Appendix  B)  and  the  function 
G™( *)  is  defined  by 


(C-l) 


(C-2) 


Interchanging  the  order  of  summation  and  integration  in  (C-l)  we  obtain 


P(e;p [space)  = \ exp 


[-  yjlSSs  [; 

n=0  m=0  L 


where 


= l expfc]  4 ^ 


1 (ft) 


(C-3) 


(C-4) 
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Using  (C-2)  in  (C-4)  and  interchanging  the  order  of  summation  and 
integration, 

U 


m. 


/h+m-jN 

1 r i 

2y  1 

\ m ) 

j! 

4(>-p2)J 

J y exp 

— oo 

Oj(l-p)J 

(C-5) 


Making  the  substitution  x=-y  in  (C-5)  and  using  [9,  eq.  3.381.4]  we  can 
evaluate  the  integral  in  (C-5)  to  obtain 

V ■ (¥)  4Z(Tj)  (•&)'  • <C-6) 

j=0 

Equation  (C-6)  can  be  summed  using  [12.P.17]  ar.d  [lO.eq.  15.4.1]  to  yield 

V ■ (¥)  4 C)  2^1  (’n*^  T+p)  . (c-7) 

Substituting  (C-7)  into  (C-3)  and  rearranging  terms  gives  (26a)  of  the 
main  text. 

C.2  EVALUATION  OF  (23b) 

Using  (20b)  of  the  main  text  in  (23b)  of  the  main  text. 


P(e;  p|mark)  = f ^ exP 


0 °d 


hl  + h2 


1 -0 


J CD  OO 

EE^r’ 

^m,n 

&($(*?♦» 

n=0  m= 

n 

exp 

■0 

Gn 

m 

4y 

dy. 

.Oj(l+p). 

2n  2% 
.°d(i-p  )J 

(C-8) 


Interchanging  the  order  of  summation  and  integration  in  (C-8)  we  obtain 

.2,  .2- 


P(e; p|mark)  . V«P  [-  ^]EEs iSf  [f  (ft)  (»?♦'!)] 


m 


n=0  m=0 


(C-9) 
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where 


m 


(C-10) 


(c-11) 


Using  (C-2)  in  (C-10)  and  interchanging  the  order  of  summation  and 
integration, 

m ao 

V.  ■£  hrJ  37  h-VT  / yJ  eJT-M 
j-o  L°d(1-p  )J  o L°d(1+p)J 

Using  [9,  eq.  3.381.4]  to  evaluate  the  integral  in  (C-11)  we  obtain 

m 

(c.«, 

j=o 

Equation  (C-12)  can  be  summed  in  the  same  way  that  (C-6)  was  summed  to 
yield 

Tm.n  '(t6)  °d  (T)  2F1  (-">.l;-i»-n;T^-)  . (C-13) 

Substituting  (C-13)  into  (C-9)  and  rearranging  terms  gives  (26b)  of  the 
main  text. 
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APPENDIX  D 

ALTERNATE  DERIVATION  OF  ERROR  PROBABILITY  EQUATIONS 

This  appendix  presents  an  alternate  derivation  of  the  conditional 
error  probability  expressions.  For  brevity,  only  the  case  of  a "space" 
being  transmitted  is  considered. 

The  probability  density  function  of  the  decision  variable  conditioned 
on  the  phase  difference  between  the  two  inputs  to  the  multiplier  is  given 
by  (14)  of  the  main  text. 

The  phase  difference  9 is  a random  variable  with  some  density  func- 
tion f^(B).  Thus  the  unconditional  p.d.f.  of  the  decision  variable  is 
determined  from 


p(y)  = y*  fy(y;p|9=B)f^(B)dB. 


(D-l ) 


Our  task,  then,  is  to  find  f.(e)  so  that  (14)  and  (D-l)  may  be  used  to 
find  the  pdf  of  y,  which  may  then  be  integrated  to  determine  the  error 
rate  performance  of  the  DPSK  system. 

PROBABILITY  DENSITY  FUNCTION  OF  PHASE  DIFFERENCE  9 
The  probability  density  function  of  9 can  be  obtained  from  the 
p.d.f.'s  of  . i'=l  ,2: 


f^(B)  = J"  * q (B+a)f^  (a)dc 


(D-2) 


In  reality,  the  functions  f and  f depend  upon  the  transmitted  signal 

*1  *2 

phases,  e1  and  respectively,  and  hence  should  be  written  as  conditional 


hi 


i 
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densities  and  Then  f+(B)  becomes  f, (fi|e^ , ) and, 

as  expected,  the  result  must  be  averaged  over  the  a priori  symbol 
probabilities. 

The  phases  4^  and  <j>2  are  identically  distributed  random  variables, 
with  density  function  [2],  [3] 


OP 

Va,ei}  = bkcos[k(a-e.)],  i=l,2. 


10-0^  1 5 11 


where 


6k  = ckTrr(ib  OjFjdi  k+  li  -R l) 

2 

R = uplink  signal  to  noise  ratio  (at  input  to 
bandpass  limiter) 


k=0 

k>0 


ek=  {2! 

and  jFjCajbjz)  is  the  confluent  hypergeometric  function. 

The  p.d.f.  of  the  decision  variable  conditioned  on  6j-e2  is 


then 


p(y)  = I ^y(y;p|^=B)f^(6|e1-e2)de 

J —00 

CD  CD 

= I I fytyjpU'B)^  (e+olej)^  (a|e2)dade  . 

. -cn  ' 00  ^ ^ 


(D-3) 


(D-4) 


The  p.d.f.  of  i=l,2,  is  given  by  (D-3).  We  first 
determine  the  limits  of  the  integration 

V6!9!-^  = J f^(e+a|e1)f^(a|e2)dQ. 


(D-5) 
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In  view  of  the  restricted  range  over  which  the  f^.'s  are  non-zero,  (0-5) 
must  be  treated  as  four  separate  cases  (see  Figure  D-l).  jf  _fi+n  4 _ 

* P Dj  ® '02*1 i 

or  e1-e2+2ir<e,  then  there  is  no  overlap  of  the  two  f#>*s  and  f^=0.  Jn  the 
second  case,  we  have  e^rc-B+ej+iKe^ir,  which  implies  that 

MB|ere2)  ^ f^1^+a|e1)f^(a)e2)da,  ej-e^Bsej-e^*. 

( D-6 ) 

In  the  third  case,  -B+ej-n^e^-HK-B+ej-Hr » which  implies  that 


= / 2 ^/b+qIoj)^  (a|e2)da,  ere?-2T<B<e,- 

e,-n-B  1 <:  1 £ i 


e. 


2 ’ 
i 

(D-7) 

In  the  fourth  case  or  ,<,^-2,.  there  js  „„  overlap  and 

f4>(Bje1-©2)  = o. 

Accordingly,  the  range  of  a and  B over  which  the  integration  (D-4) 
should  be  performed  is  as  shown  by  the  shaded  area  in  Figure  0-2 : 

re2+«,fira+' 

P(y)  = / / V(ale2)fy(y;PU=fi)f.  (B+a|e  )deda.  (D-8) 

To  calculate  the  error  probability,  consider  the  case  where  a space 
is  transmitted,  i.e.,  there  is  no  change  of  transmitted  phase  from  one  symbol 

to  the  next,  Bj-e^O.  Thus  an  error  is  made  if  the  decision  variable  y<0,  and 
we  have 


l 


P(e;p (space)  = | p(y)dy 

.0  . e?+jt  r e,-a+7i 


r 0 j -o+7i 

J % tQlej)f(y;pU=B)f  (B+aleJdBdc 
>2_u  * *1 

(D-9) 


dy 


D-3 
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c)  Case  3 d)  Casc  A 


FIGURE  D-l 


DETERMINATION  OF 
FOR  PROBABILITY 
A FUNCTION  OF  fi 


INTEGRATION  LIMITS 
DENSITY  OF  * AS 
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Putting  (14)  of  the  main  text  into  (0-9)  we  have,  after  some  algebraic 
manipulation. 


P(e;p|space)  = — exp 


-K 


CD  CD 


"!  Ul+p  )(hl+h2}_ 


m 


where 


and 


Y = 


2phjhz 

~T7 


r at  l G 

m -*y 

.°d(1-p). 

n 2,.  F 
od(l-p  ) 

JrO^+Ti  r Oj-a+n 

I j eY  C0stl|^l+z  coseJniri-Zcosp1n 

a •'n 


• ( 3+«  1 0 j )dpda 


( D- 1 0 ) 


2hjh2 

L 5 — 2 

hi+h2 


Let  us  denote  the  double  integral  in  (0-10)  by  1^.  Then 


IS  = 
j2 


re2+v  s 

= I (n|ej)I3(a)da 

*^flo  “11 


( D- 11) 


where 


/•  0 j-a+Ti 

(a)  = I expj^Ycose]  [l4Zcose]m[l-Zcose]nf^  (p+alfljJdB. 

-V-a-n  1 


'1 


We  can  use  the  binomial  theorem  to  write  (D-12)  as 
m n 


i: 


(D-12) 


(D-13) 


l=0  r=0 


0-6 
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where 


6, -a+H 


rru" 

•^0j-a“n 


Jc0SB„  M 


cos  Bf  (B+aleJdB 
*1  1 


(0-14) 


and  M = £ + r. 


Substituting  (D-3)  into  (D-14)  we  have 


.0, -a+n 


V(M)  = 2"  2^bk  f eY  COSecosMecos[k(B+a-e  )]dB.  (D_i5) 

k=0  Je}-a-v  1 

Using  [9,  1.320.5  and  1.320.7]  we  can  expand  the  power  of  the  cosine  in  terms 
of  cosines  of  multiple  arguments 


cosMc 


M-u(M) 

2 


2M  VotM)  ([M-u(M)-2p]/2  ) C0S{f2M+u(M)] c}  (D-16) 


where 


e = f 1 > ^ =0 

11  ( 2,  v>0 

»<«>  ■ {?:  5 o 


q even 
q odd 


Using  the  identity 


COS  Z COS  C = ljcos(z+C)  + cos(2-o] 


(D-17) 


in  connection  with  (D-16)  allows  us  to  reduce  the  integral  in  (D-15)  to  the 
form  [13] 


v = J eYcoSacosNa  da  = 2uIn(Y) 
▲ N 


(0-18) 


where  $ is  an  arbitrary  angle. 


D-7 
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Substitution  of  (D-16),  (D-17),  and  (D-18)  into  (D-15)  yields,  after 
some  algebraic  manipulation. 


M-u(M) 

V(M)  53  C0^MQ-e  Be  b 

* k=0  p=0  Vh,(M)  k 


[Vu(M)-k(Y>  + ^w+uCHj+k^)] 


([M-u(M)-2p)/2j 


(D-19) 


where  bk,  C|(,  and  u(k)  are  as  defined 


and  (D-3)  in  (D-ll)  yields 


Previously.  Using  (D-13),  (0-19) 


l+r-u(i+r) 


m n 2 


£EEEE  E 

p=0  k=D  1=0  r=0  P =0 


hMml  MA  »_nr7l+r 
2l+r  1 1 A r J ( 


i+r 

l+r-u(i+r)-2u 

2 


* [^uMi+rJ-k^  + ^u+uti+rJ+k^U  bpbk  f cos  ^(ft-ejjJcos^kfa-ejjJda. 

•'a 


(D-20) 


The  integral  in  (D-20)  is  equal  to  2*  for  p=k=0,  , for  p=k*0,  and  disappears 


for  p^k. 


D-8 
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Thus  we  have 


£+r-u( t+r) 
m n 2 


4 E E E E i(X)2ltr(-Dr 

^ * k=n  f=n  r=n  .,=o  ti,vlAr'  0i+r 


k=0  £=0  r=0  u=0  cj( 


^(i+r“u( £+r)-2g]/2  )|j  2p+u( £+r)-k^  4 1 2u+o( t + r)  + k^ . (D-21) 


Making  use  of  the  property  of  the  binomial  coefficients. 


k J = 0 if  b<0  or  b>a , 

to  extend  the  limits  of  the  finite  sums  to  infinity  and  then  applying  Bailey' 
theorem  [11, pp.  58-59]  to  the  sums  over  £ and  r,  we  can  write  (D-21)  as 


(-D5 


[W)-k(v)tW)n<Y>]  i (tXs-t)  <-1)t- 


(D-22) 


It  can  be  shown  [12,  p.  17]  that  the  summation  over  t in  (D-22)  can  be 
expressed  as  a Gauss  hypergeometric  function  with  a negative  numerator 
parameter.  Making  use  of  this,  we  have 


s-u(s) 
n 2 


Is  - 1 E T T I 5 Vz  f b>V»i!)  , ,)»(s) 

■[I2u*u(s)-k(V>*I2l.-tu(s)+k(’,)] 

*(s)2Fl^m,'S;  n's+1; 


(D-23) 


D-9 
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where  we  have  again  used  the  properties  of  the  binomial  coefficients  to 
introduce  finite  summation  limits  and  have  used  (-1)S  = (-l)u^. 


Letting  1^  denote  the  first  integral  in  (D-10),  from  the  definition 


of  the  Dolynomial  G™(x)  we  have 


/ 

,S  = £ /n+m-j 


V m 

j=G  X 


)* 


tK- l)j 


4 

T° 

L^J 

L p 

yjdy.  (D-24) 


The  integral  in  (D-24)  may  be  evaluated  by  making  the  substitution  x = -y  and 
using  (9,  eq.  3.381.4]  to  obtain 


(D-25) 


Using  [12],  we  can  sum  (D-25)  to  yield 


(0-26) 


Substituting  (D-23)  and  (D-26)  into  (D-10),  and  rearranging  terms  for 
clarity,  yields: 


I • i 


i 


D-10 


I 
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P(e;p|space)  exp  - 2- 


s-»(s) 


•p  n=0  m=0  k=0  s=0  p= 


ckcp+u(s)^  ^ ^ 


(?) 


(s)(ls-U(L)-2u]/2))(f)(~2)  [^)Rd  ('+*2)] 


. to  R2(1+x2)l  ! |2pXRd  ) + ! ^2pXRd 

|^4  \1  -P / Rd(1  }J  VofsJ-k  \^?J  I2ptu(s)+k^^' 

* 2Fl  ^-n.li-n-m;  2Fl^“m,"s;  s+1^"1) 

* [lFl(l;  k+lj-Ru)]2 


(D-27) 


where 


cq  "(*!  q>0 

u(q)  ={°; 


q even 
q odd 


and  where  we  have  defined  the  signal-to-noise  ratio  parameters 
2 2 

R^  = h.|  = direct  channel  SNR  into  the  multiplier 


2 "2 

A = — £■  * power  imbalance  between  multiplier  inputs 
hl 

and  we  have  used  the  relation  (2z)!  = r(2z  +1)  = - r^  + |jr(z+l)  to 
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APPENDIX  E 

ALTERNATE  FORMS  FOR  ERROR  PROBABILITY  EQUATION 


The  error  probability  equations  (26a)  and  (26b)  of  the  main  text 
may  be  written  in  several  alternate  forms.  Although  these  alternate 
forms  do  not  appear  to  be  any  more  computationally  efficient,  they  are 
presented  here  for  the  benefit  of  those  readers  who  may  wish  to  pursue 
further  the  mathematical  properties  of  the  error  behavior  of  a DPSK 
system.  For  brevity  we  present  only  the  forms  fora  "space"  being 
transmitted;  the  forms  for  “mark"  can  easily  be  written  by  analogy. 

First,  we  consider  the  inner-most  summation  of  (26a): 

s-u(s) 

2 

A = Cv+v)(s)([s-v)(s)-2p]/2)  [I2y-h)(s)-k(Y)  + I2M4u(s)+k(Y)]  (E_1) 


2 2 

where  Y = 2pXR^/( 1 -p  ).  The  expression  in  (E-l)  may  be  summed  as  shown 


below. 


From  [10,  eq.  9.6.29]  we  find  that 

i K,-p<2>  + (?)  'v.p*<2>  ♦(?)  * 


* Iv+p(z )j»  P 0 » 2,... 


where 


(E-2) 


I . (z)  - 


,P  v 


I (z)  . 


For  p=0,  l*u,(z)  = !v(z).  Also,  from  [10,  eq.  9.6.6] 


]-n(l)  = !n(z) 


(E-3) 
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Using  (E-7)  in  (E-l),  we  find  that 


A = 2s+1l[s)(Y)  . 

Putting  (E-8)  into  (26a)  of  the  main  text,  we  have  the  alternate 
expression 


(E-8) 


P(e;p  ,x|  space)  = n(^)  exp 


\ 

2Rd 

) exp 

" , 2 

1-P 

“ • ” n , ..u(s) 

n=0  m=0  k=0  s=0  n,!n![I\~r/J 


2 

• 2E1(-n».-s;n-s+1;-1)|1Fi(|-;k+l;-R^j 

for  the  error  probability  given  that  a "space"  was  transmitted. 

Another  formulation  for  the  error  probability  may  be  written  by 
using  the  relationship  [14] 


(E-9) 


r . i 2v  exp(-fr/2)r(n+l )2 

1F1(n+i;2n+2;-R^)  = — u 


2n 


w 


['ir)*  '„«(?)]•  (E-'°) 


Setting  k/2  = n+£,  the  confluent  hypergeometric  function  in  (26a)  of 
the  main  text  may  be  expressed  as 

r ( p2j  _ exP(-Ru/2)r(JTI)2k'1r  Au\] 

1F1(  2*k  Ru)  - 7R2j(k-l)/2  l^tk-l )/2\T/+  I(k+1)/2VT/J* 

(E-ll) 

Putting  (E-ll)  into  (26a)  yields  a second  alternate  formulation  for  the 
error  probability 


E-3 
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P(e;p,x| space)  = irR^  (-y|-)  exp 


2R 


1-P 


exp 


K] 


Squill 


Ck  Cyhj(s) 


(-D 


u(s) 


min! 


n=0  m=0  k=0  s=0  y=0 


(r)(-)(^ 


• > fek2  »«2>r  [i  <’«2>]n 


r2pxR 


2p+u(s)+k\  2 

l-p 


2x1 


(k+l)/2 


m 


o 

2F1(-n,l  rn-m;  ) ^Fjt-m.-sjn-s+l  ;-l)  . 


(E-12) 


A third  alternative  form  may  be  derived  by  using  both  (E-8)  and 
(E-ll)  in  (26a)  to  yield 


T 2R^  "I 

P(e;p ,x|space)  = nR*  exP M exp["Ru] 

1 ”P  J • 


CD  CO  f") 


• EEEEth  (n;) (") 

n=0  m=0  k=0  s=0  u A ' 


• [i  (i^Rd  ht  & (w)  Rd  (i+i2>]n 


['ik-l)^"?1)*  1 ( k+1  )/2 ) 
2 


• 2Fi(“n»1  yf^-)  2Fi(-m,-S;n-s+l;-l) 

where  we  also  have  used  (— 1 )s  = (-l)u^. 


(E-13) 


i ■ rii'  tmriiMi  ... 


K- 

& 
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An  interesting  fourth  transformation  arises  from  the  use  of  Bailey's 
theorem  [11]  to  replace  the  infinite  sums  over  n and  m in  (E-13)  with 
one  infinite  sum  and  one  finite  sum.  Letting  v=n  and  w=n+m  we  obtain  from 
(E-13),  after  a little  algebra. 


P(e;p,x|space)  = nR*  exp["Ru]  exp [' 


w 


k=0  w=0  v=0  s=0 

,2 


[i  tew  <~'.r 


. (J«.f  ,w(aA)  [l„. „„(!?).  I, /f)f 


♦ 2Fj(-v,1;-w;  yf^-)  ^(v-w.-s  ;v-s+l  ;-l) . (E-14) 

The  form  in  (E-14)  begins  to  resemble  summation  formulas  for  generalized 
hypergeometric  functions,  for  example  [15]  or  [16].  This  is  an  area  of 
mathematics  which  is  not  yet  fully  explored;  hence  the  implications  of 
forms  such  as  (E-14)  cannot  be  fully  stated  at  this  time. 

Other  forms  can  be  written  by  applying  various  transformation 
formulas  to  the  Gauss  hypergeometric  functions  in  (26a),  (26b),  (E-9), 
(E-12),  (E-13),  and  (E-14).  Examples  of  applicable  transformation  can 
be  found,  for  example,  in  chapter  15  of  [10].  It  seems  of  little  value, 
though,  to  write  out  in  full  the  multitude  of  forms  thus  derivable. 


E-5 
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APPENDIX  F 

A FURTHER  CONSIDERATION  ON  POWER  IMBALANCE 

The  power  imbalance  between  adjacent  pulses  arises  due  to  such  effects 
as  delay  line  attenuation,  phase  error,  and  intersymbol  interference. 

The  power  imbalance  may  also  arise  from  sources  other  than  the  receiver 
itself  such  as  fading.  The  effect  of  delay  line  attenuation  or  fading 
on  power  imbalance  is  straightforward  and  need  not  be  mentioned.  The 
phase  error  can  arise,  for  example,  from  an  improper  delay  length 
in  the  delay  circuit  of  a phase  detector.  The  intersymbol  interference 
is,  however,  a complicated  problem  resulting  from  the  combined  effects 
of  the  delay  circuit  and  the  transfer  function  (filter  characteristics) 
of  the  system  under  consideration. 

In  this  appendix  we  want  to  enhance  our  understanding  of  the 
possible  relationships  between  power  imbalance,  intersymbol  interference, 
and  phase  error  associated  with  non-ideal  delay  lines  in  the  differential 
phase  detector. 

Hubbard  [17]  analyzed  the  effect  of  intersymbol  interference  on 
the  probability  of  error  (for  the  case  of  no  noise  correlation)  under 
the  assumption  that  the  intersymbol  interference  comes  only  from  adjacent 
pulses,  and  showed  that  the  probability  of  error  is  a function  of  power 
imbalance  (between direct  and  delayed  channel)  caused  by  intersymbol 
interference  [17,  eq.  14]. 

Another  way  of  viewing  the  power  imbalance  is  to  relate  it  to 
phase  error.  The  intersymbol  interference  is  assumed  to  manifest 
itself  in  the  form  of  a perturbation  in  the  phase  of  the  signal  at  the 
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sampling  instant.  More  precisely,  the  intersymbol  interference 

introduces  a phase  shift  A<#>  so  that  the  phase  change  in  a time  slot  is 

<j>±A$  instead  of  <p  U=0  for  a "space"  and  $=u  for  a "mark").  The  value  of 

depends  on  the  details  of  the  signal  waveform  and  the  filter 

characteristics  of  the  system.  The  determination  of  the  A<f>'s  (for  a 

particular  system)  which  describe  the  intersymbol  interference  phenomena 

in  a complicated  manner  is  a rather  difficult  problem,  and  even  if  it 

is  theoretically  possible  (see  [18]  for  example),  the  actual  measurement 

of  this  quantity  will  be  a formidable  task.  For  this  reason  we  will 

establish  an  analytical  background  to  replace  A<t>  by  an  easily  measurable 

2 

quantity,  power  (or  equivalently  SNR)  imbalance  X which  is  used  as  a 

2 

basic  parameter  throughout  the  report.  By  relating  A4>  to  X , all  the 

possible  degradation  factors  mentioned  above  (attenuation,  fading,  phase 

error,  and  intersymbol  interference)  which  otherwise  appear  to  be 

2 

different  attributes  are  merged  into  one  quantity  X . The  statistics 
2 

of  X may  then  be  experimentally  determined  for  a particular  system. 

2 

We  now  show  the  relation  of  A$  to  X . To  do  this  we  note  from 

the  definitions  following  (14)  in  the  main  text  that  the  parameters  h^ 

3 

and  h^  are  all  that  we  need  in  describing  the  error  performance  in  terms 
of  power  imbalance.  For  simplicity,  consider  the  case  of  no  noise 
correlation  (p=0).  We  may  write  then 


HP 

2h3- 

hi  + h2  * 2hjh2 

cos$ 

( F-la) 

2h4  = 

hl  + h2  ' 2hlh2 

COS<J> . 

( F-lb) 

F-2 


Let  us  define  the  average  SNR  (or  equivalently  average  power)*  by** 


h*  + h* 

K - * 


and  a differential  phase  (or  phase  error)  A$  such  that 


COSA<f>  = 


2hjh? 


2 2 
hl+h2 


Then  (F-l)  can  be  written 


? 7 7 7 2h.  h» 

hj  + h^  + (hj+hp  -2 1 2 cos«j> 

hl+h2 


? 2 

2R  + 2R  cosA<f>  cos$ 


2 2 2 
= 2R  - 2R  cosA$  oos<J>. 


Also,  defining  the  power  imbalance  x by 


2 A 2 

x = 7 

hl 


we  can  write  (F-3)  as 


COSA$  = 


2 _ 1 -s i nA4> 
l+sinA<f> 


(F-2) 


(F-3) 


(F-4a) 


( F-4b) 


(F-5) 


*In  binary  communication  systems,  error  performance  of  a correlation 
receiver  (optimum)  depends  only  on  the  average  signal  energy  E=(E0+Ej)/2 


where  EQ  and  Ej  are  energy  associated  with  bits  "0"  and  "1"  respectively 
provided  the  level  of  noise  spectral  density  is  the  same  in  two  inputs 
and  the  two  input  signals  are  uncorrelated  (19,  p.  163]. 

**  This  definition  of  R differs  from  that  used  in  Appendix  H of  this 
report. 


■ f 
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We  note  that  (F-l)  describes  H3  and  in  terms  of  hj  and  h2  (two 
different  input  SNR's),  while  (F-4)  describes  the  very  same  parameters 
by  R and  (constant  average  power  and  differential  phase  error).  We 
may  call  the  former  a "power  imbalance  model"  (since  h^h^  in  general), 
and  the  latter  a "phase  imbalance  model."  The  important  fact  is  that 
the  two  models  are  equivalent  under  the  constraints  of  (F-2)  and  (F-5). 
A geometrical  interpretation  will  be  more  instructive  in  grasping 
physical  understanding  of  the  relationships  between  the  parameters 
involved. 


Let  us  consider,  for  simplicity,  the  case  that  a "mark"  is  trans- 
mitted; <J>=ei_62=1T*  Then  since  cos4>=-l  and  sin4>=0,  we  may  write  (F-l) 
and  (F-4)  as 


h3  . hj-hj 

H,  = h,*h, 
4 12 


and 


H-  = yR2+R2-2R-R  cosA<|> 


(F-6a) 
( F-6b) 

( F-7a) 


H4  =^R2+R2+2R*R  cosA4>  . (F-7b) 

We  now  show  that  the  geometry  in  Figure  F-l  satisfies  all  the 
relationships  given  by  (F-6),  (F-7),  and  also  (F-2)  and  (F-5).  Noting 
that  CP=EP=hj,  we  have  (F-6a),  while  (F-6b),  (F-7a)  and  (F-7b)*  are  obvious 
from  the  figure.  Equation  (F-2)  follows  from  the  fact  that  CD2=CP2+DP2=  h2+h2; 
and  CfD2+AD2=2CD2=AC2=(2R)2.  The  validity  of  (F-5)  can  be  shown  as  follows: 

X = = — = tangent  of  angle  DCP  = tan  (45°-a<j>/2)  . 

nl  CP 


1 


: 

\ 


* (F-7a)  and  (T-7b)  are  from  the  law  of  cosines. 


F-4 


6a  = oc  = ob=  r 
ab=be=h3 
bc  = h4 

CP  = h, 

BP=DP  = h2 

FIGURE  F-l  GEOMETRICAL  REPRESENTATION  OF  PARAMETERS  H3  AND 

IN  TERMS  OF  hj  AND  h£  (POWER  IMBALANCE  MODEL),  AND 
R AND  (PHASE  IMBALANCE  MODEL) 
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Thus 

X2-  tan2{45°-A4/2)  - , l^°s(90°-af) 

cos^(45°-a*/2)  l+cos(  90°-A<f>) 

1 -si  nA4> 

1 +s  i nA4> 


( F-8) 


Equation  (F-8)  or  (F-5)  is  significant  in  that  it  bridges,  under  the 
constraints  of  constant  average  power  (F-2),  the  gap  between  the  power 
imbalance  model  and  the  phase  imbalance  model.  The  relationship 
between  X and  A#  as  a function  of  the  degree  of  imbalance  can  also  be 
shown  from  the  figure.  When  point  "P"  is  close  to  point  "0",  we  can 
see  that  hj=h£=R  and  A<J>=0  which  represents  an  ideal  situation.  As  the 
point  "P"  moves  toward  point  "D"  along  the  arc  OPD,  the  value  of  X=  h^/hj 
decreases  (power  imbalance  increases)  and  accordingly  A$  increases. 

As  "P"  approaches  "D",  X becomes  zero  (large  power  imbalance)  and  A$ 
becomes  90°.  Thus  we  can  see  that  change  of  X from  1 to  0 corresponds 
to  change  of  A<J>  from  0°  to  90°. 


j 
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APPENDIX  H 

ERROR  BEHAVIOR  OF  A BINARY  DPSK  SYSTEM 
DUE  TO  INTERSYMBOL  INTERFERENCE  AND 
CORRELATED  NOISE* 

It  is  well  known  that  the  classical  result  of  the  error  probability 
for  differentially  coherent  detection  of  binary  PSK  applies  only  to 
the  ideal  situation  where  the  received  symbol  signal  energy  from  pulse 
to  pulse  is  assumed  equal  and  the  noises  at  the  sampling  instants 
uncorrelated.  The  purpose  of  this  appendix  is  to  show  the  error  behavior 
of  a differential  phase  shift-keying  (DPSK)  system  under  practical 
assumptions  where  the  signal  powers  at  sampling  instants  between  the 
adjacent  pulses  are  unequal  and  the  .noises  cbrrelated.  Error  performance 
is  presented  graphically  for  different  levels  of  power  imbalance.  A 
significant  result  observed  is  that  the  probability  of  error  is  indepen- 
dent of  noise  correlation  for  all  degrees  of  intersymbol  interference 
(power  imbalance).  This  result  (based  on  our  computations)  is  clearly 
a departure  from  the  previous  beliefs. 

ANALYSIS  MODEL 

The  binary  DPSK  system  under  consideration  is  depicted  in  Figure 
H-l.  The  primary  objectives  of  this  appendix  are  the  considerations  of 
the  error  performance  calculations  under  the  assumptions  where  (SNR) ^ (SNR)^ 
and  the  noises  n^(t)  and  n2(t)  are  statistically  dependent  (see  Figure  H-l). 
The  situation  where  (SNR)j  is  different  from  (SNR)2  is  brought  about  by 
unequal  signal  powers  at  the  phase  detector  of  the  DPSK  demodulator,  due  to 

* This  appendix  is  based  largely  on  [5]  with  a change  in  the  definition 
of  the  SNR  parameter.  Minor  typographical  errors  in  [5]  have  also  been 
corrected. 
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such  phenomena  as  intersymbol  interference,  phase  errors  in  the  delay  circuit 
and/or  delay  circuit  attenuation.  Thus,  the  error  performance  analysis  based 
on  the  unequal  signal  powers  over  two  consecutive  symbol  signal  pulses  would 
constitute  an  upper  bound  of  error  probability  in  the  presence  of  intersymbol 
interference.  The  noise  correlation  is  a practical  assumption  since  the 
noise  is  necessarily  bandlimited  in  practice. 

In  Figure  H-l,  the  two  inputs  at  the  phase  detector  (multiplier), 
uj(t)  and  u2(t),  represent  the  noisy  received  waveforms  for  two  consecutive 
source  symbols: 


Uj(t)  = SjU)  + nx(t) 

(H-l) 

u2(t)  = S2(t)  + n2(t) 

(H-2) 

where  S^(t),  i=l,2,  are  the  bandpass  information  carrying  signals;  and 
n-(t),  i=l,2,  bandpass  noises.  Assume  that  the  "present"  and  the  "preceding" 
source  symbols  are  identified  with  carrier  phases  Bj  and  e^,  respectively. 
Then  we  may  write 

3j(t)  = /2P^cos  (iot  - Bj)  (H-3) 

S2(t)  = ^2P^cos  (w(t-T)  - e2] 

= /2P2  cos  (wt  - e2)  (H-4) 

where  the  carrier  frequency  is  assumed  to  be  chosen  such  that  wT  = 2nk, 
k integer,  and  Pj  and  P2  are  the  carrier  powers  at  the  phase  detector.  We 
shall  be  primarily  concerned  with  the  case  where  Pj  f P2. 

The  bandpass  Gaussian  noises  are  expressed  in  the  forms: 

nl ( t)=Xj (t)  cos  lot  + Yj  (t)  sin  iot  (H-5) 

n2(t)  = Xj(t-T)  cos  [to(t-T))  + Yj(t-T)  sin  [io(t-T)] 

= X2(t)coswt  + Y2(t)sinut  (H-6) 
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and  the  noise  correlation  is  defined  by 

E(n1(t)n2(t) } = EfnjUJnjU-T)}  = o2p(T)  (H-7) 

where  p(T)  is  the  normalized  correlation  with  T equal  to  symbol  duration, 
and  the  noise  power  is  given  by 

o2  s E(n2(t) } = E{X2(t) } = E { Y?( t) } ; i=l,2.  (H-8) 

The  decision  variable  Y(t)  at  time  t is  the  output  of  the  zonally 
low  pass  filtered  version  of 
X(t)  = uj(t)  x u2(t) 

-[VH^cos  (wt  - 6j)  + Xj(t)  cosut  + Yj ( t)  sinut] 

. [/2P_2  cos  (wt  - e^)  + X2(t)coso)t  + Y2(t)  sinut]  . (H-9) 


Since  6^  - e2  is  either  0 or  +ir  in  a binary  DPSK  system,  the 
binary  decision  is  based  on  the  comparison  of  the  decision  variable  Y(t) 
with  "zero"  threshold,  for,  when  noises  are  assumed  to  be  absent  at  the 
input,  the  decision  variable  is  given  by 

Y(t)  * 4^cos  (ej  - e2)  . (H-10) 

DECISION  VARIABLE  STATISTICS 

The  problem  of  obtaining  the  probability  density  function  (pdf) 
of  the  lowpass  filter  output  for  the  system  model  that  fits  our  situation, 
depicted  in  Figure  H-l,  was  solved  by  Miller  and  Lee  [4]  in  great  generality. 
Our  need  here  is  a special  case  of  the  problem  treated  in  [4]  . From 
14,  eq.(25))  we  have  the  pdf  for  the  decision  variables  y(t)  as  follows: 
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where 


H3  ^ 2(l+p)  [hl  + h2  + 2hih2cos  ^ei  " e2j] 
H4  ” 2(  1-p)  [hl  + h2  " 2hih2cos  ^ei  ' v] 


g:(2)  a 5 


( m+n-k\ 

h\  n /k! 


and  where 


? P. 

hj  ^ -js  (SNR).;  i=l ,2 

o 


(H-12a) 

(H-12b) 
( H-l 2c ) 

( H -1 2d ) 


and  P. , a , p=p(T),  are  the  signal  power,  noise  power  and  noise  correlation, 
respectively,  as  defined  earlier. 

If  we  assume  that  the  reference  phase  e2  is  0,  then  6j  is  either 
0 or  n,  depending  upon  whether  the  symbol  following  the  reference  symbol  is 
0 (space)  or  1 (mark),  so  that  = 0 or  6j  - = it.  The  conditional 

pdf's  are  thus  obtained  as  follows: 

f(y;P|space)  = f (y;P  | ej-e^O)  . (H-13a) 

f(y;p|mark)  = f (y;p  | ej-e^ii)  . (H-13b) 


CONDITIONAL  ERROR  PROBABILITIES 

From  (H-13a)  and  (H-13b)  one  can  obtain  the  conditional  probabilities 
of  error  from  the  expressions  q 

P(e;p|space)  = Prob{y<0| ej-e2=0)  = j f (y ;p | ej-e2=0)dy  (H-14a) 

P(e;p|mark)  = ProWy>0|e1-e2=7r)  =J  f(y;p|e1-e2=7i)dy.  (H-14b). 


and 


Carrying  out  the  integrations  indicated  by  (H-14a)  and  (H-14b)  using  the 
density  functions  given  in  (H-ll),  one  obtains  the  following  results: 
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and 


P(e;p|space)  = ^-expj^ ^r(hj  + h2  ' 2Phjh2)j 

V*  V V'2'^m+n)  / n+m-k  ^ 

h£oikmlnl  \ * / 


. [K&K  ^>T[K&) <h:  - vfW  tH-15., 


P(e;p|mark)  = exp  J ^ (hj  + h2  + 2ph^h2) 


hi 

- - ^2-(^n)  / n.k\ 

n=0  m=0  k=0  m'ni  \ n / 

<-»»> 


It  is  interesting  to  observe  a "symmetry"  property  of 
P(e;p|space)  = P(e;-p|mark)  . 


(H-16) 


ERROR  PROBABILITY  EXPRESSIONS 

To  compute  the  error  rates  using  the  equations  (H-15a)  and  (H-15b),  we 
need  to  define  an  important  variable.  In  a DPSK  system,  a single  symbol 
decision  is  made  using  two  symbols.  Since  we  are  considering  a situation 
where  each  symbol  signal  energy  is  not  equal  (due  to  intersymbol  interference, 
for  example),  it  is  appropriate  to  define  the  signal-to-noise  power  ratio 
per  pulse  (or  symbol)  as  follows: 

= hj  = direct  channel  SNR.  (H-17) 

. ■■  (a  classical  case),  there  is  no  difference  between  this  defini- 

• the  conventional  SNR  definition  per  symbol. 
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Define,  further,  an  SNR  difference  measure  by 
2 (SNR), 

X T5HKJT 


(H-18) 


2 2 

In  terms  of  R and  A , the  conditional  error  probabilities  of  (H-15a) 


and  (H-15b)  are  given  by 


P(e;p;A|space)  = 


' ^ exp[' ' 2ol)j 

•SSS 


(H-19a) 


P(e;p;A|mark)  = ^ exp  + x2+  2Px 


>-(m+n) 


n=0  m=0  k=0 


?(T‘)ter 


[r  "»’f  [i  "-’tfcj 


(R  ) . (H-19b) 


The  total  unconditional  probability  of  error  is  the  weighted  sum  of 
(H-19a)  and  ( H— 1 9b ) given  by 

P(e;p;X)  = PsP(e;p;x|space) + PMP(e;p;X|mark)  (H-20) 

where  P$  and  PM  are  the  aj)riori  probabilities  of  space  (0)  and  mark  (1). 


r 


J.  S.  LEE  ASSOCIATES,  INC. 


ERROR  BEHAVIOR  UNDER  SPECIAL  CONDITIONS 


If  x^=l  is  substituted  into(H-19a)  and  (H-19b),  we  obtain 


P(e;p |space)  = ^ exp(-R^) 

P(e;p|mark)  = exp(-R^) 

and  the  unconditional  error  probability  of(H-20)  becomes 

2 

P(e;p)  = |[l  + (PM  - Ps)p]e'R. 


(H-21a) 

(H-21b) 


(H-22) 


This  is  the  identical  result  given  by  Lee  and  Miller  [7].  The  significance 

of  the  result  given  by  (H-22)  is  that  the  probability  of  error  depends  on 

noise  correlation  p only  if  a priori  probabilities  are  unequal.  When 

Pc  = Pu  = ? . (H-22)  reduces  to  the  classical  error  rate  expression  of 
b n c 

p(.)  - kr2. 

Now,  let  us  assume  that  mark  and  space  are  equi -probable: 

P^  = P^  = When  one  observes  the  conditional  probabilities  of  error  as 
given  in  (H-19a)  and  (H-19b),  it  appears  certain  that  the  total  unconditional 
probability  of  error  given  in  (H-20)  still  depends  upon  noise  correlation  p. 
In  fact,  it  has  been  remarked  in  the  previous  publications  [6],  [7]  that 
the  error  probability  of  a binary  DPSK  system  depends  on  the  correlation 
if  there  is  intersymbol  interference  and  that,  in  the  absence  of  inter- 
symbol interference, the  error  probability  is  independent  of  the  noise  cor- 
relation provided  that  P^  = P^  = 

We  have  computed  the  probability  of  error  expression  (H-20)  with 
1 2 

PM  = ^ for  various  values  of  A (different  degrees  of  intersymbol 

interference)  and  p.  We  have  observed  the  surprising  results  that  the  com- 
puted probabilities  are  independent  of  noise  correlation  p for  all  values 
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2 

of  A considered!  The  conditional  probabilities,  however,  were  dependent 

upon  noise  correlations.  When  these  conditional  error  probabilities  were 

added  with  equal  weighting  (PM  = = resu^t  converged  to  the 

values  of  error  probability  for  the  case  p=0.  To  show  the  error  mechanism, 

we  have  plotted  in  Figures  H-2  and  H-3  some  specific  cases  of  conditional 

error  probabilities.  Figure  H-4  shows  the  error  probabilities  of  a 

DPSK  system  for  A = 0 dB  (no  intersymbol  interference),  -1  dB,  -2  dB,  and 
2 

-3  dB.  Note  that  A = 0 dB  corresponds  to  the  ideal  classical  case.  Our 
results  indicate  that  the  probability  of  error  for  a binary  DPSK  system 
depends  on  noise  correlation  p only  when  the  message  symbol  probabilities 
are  unequal.  Whenever  the  prior  probabilities  are  equal  the  error  rate 
is  independent  of  noise  correlation  whether  or  not  there  is  intersymbol 
interference.  As  stated  earlier,  this  is  a departure  from  the  previous 
beliefs.  It  must  be  stressed,  however,  our  findings  are  based  on  the 
computational  results.  The  mathematical  complexities  of  the  error  rate 
expressions  did  not  lend  themselves  to  an  analytical  verification  of 
the  computational  results  observed. 


ALTERNATE  FORMS  FOR  THE  ERROR  RATE  EXPRESSIONS 


The  error  rate  expression  in  ( H-l 9a ) and  (H-l 9b)  may  be  written  in 
several  alternate  forms.  First,  we  can  use  the  relation 


(nr) 6 *>» 

to  replace  the  finite  upper  limit  of  the  summation  over  k in  ( H-l 9a ) 
and  (H-l 9b)  by  infinity.  If  we  then  let  £=  n-k  in  (H-l 9a ) we  obtain 


(H-23) 


,-mi  V ' - iV- 
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P(e;p ,x| space)  = -^c-  exp  — ( 1 +A 2-  2pA) 

1-p 


m=0  k=0  1=0 


(H-24) 


A similar  form  can  also  be  obtained  from  (H-19b)  for  the  case  of  a 


mark  being  transmitted;  for  brevity  we  will  not  go  into  the  details  of 
that  case  in  the  following  discussions. 


Interchanging  the  order  of  summations  over  k and  i.  in  (H-24), 
the- summation  over  k is  recognized  as  a confluent  hypergeometric 


function: 


£■  A a 


(H-25a) 


Applying  Kummer's  transformation  [10,  eq.  13.1.27]  to  (H-25a),  we  obtain 

r>  “■>[t  (H-25b> 


We  next  interchange  the  order  of  the  summations  over  t and  m and 
recognize  the  summation  over  m as  a confluent  hypergeometric  function: 


= iFi[£+i*i;  R4  ] - 


(H-26a) 


Applying  Kummer's  transformation  [10,  13.1.27]  to  (H-26a)  we  find  that 
(H-26a)  can  also  be  written  as  a Laguerre  polynomial  [10,  13,6.9]: 


exp  [V(i^)(1+X)2]  Li["  t (r?p){1+x)2J 


( H-26b) 
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Using  (H-25b)  and  ( H— 26b)  in  (H-24)  we  find  the  alternate  form 
P(e;„.x|space)  - «p[-  £ (l^) »->■)'] \ [-T^)'1",2] 


j Fj  |t ;«.+!; 


R2  (1-X)2 
T l-p 


(H-27a) 


Similarly  for  mark  we  can  derive  the  form 


P(e;p  ,A|mark)  ■ ^ axp  [-  £ (Hx)2]^  ^ (^tl-x)2?  (|^)(1*x)2] 

A=0  L J J 

• lFl[t;l+1-  - T %^2]  • (H-27b) 

The  forms  (H-27a)  and  (H-27b)  are  more  computationally  efficient  than 
(H-19a)  and  (H-19b)  and  are  the  forms  implemented  in  the  computer  program 
listing  contained  in  Appendix  I. 

Returning  to  (H-24),  we  can  write  another  alternate  form  for  the 
error  rate  expressions.  Using  the  relation  [2,  eq.  A. 1.44c] 


n!  = (1). 


(H-28) 


where  (z)n  is  Pochhammer's  symbol  [10,  eq.  6.1.22]  the  triple  summation 
in  (H-24)  can  be  recognized  as  a confluent  hypergeometric  function  of  three 
variables  [16,  eq.  2.13]  and  we  find  that 

P(e;p,xJ  space)  = ^ exp  - -L  (1+X2-  2pX)  ,1 ,1  ;1 ,1 ,1  ;XC  ,Y  ,Z  ) 

l-p  J s s s 

(H-29a) 
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where 


(in)2. 


Y = 
s 


r£  A-p\ 

4 \l+p) 

2d4 

* 


and 


Z 

s 4 


\A-p/ 


Similarly  we  can  derive  the  other  case 


P(e;p,x|mark)  = exp£-  (l+X2+2Px)j  3«^3)(1 ,1 ,1;1 ,1 ,1  jX^Y^) 


(H-29b) 


where 


The  forms  (H-29a)  and  (H-29b)  are  more  compact  than  other  notations, 
but  their  full  analytical  importance  is  unknown  at  this  time  since  the 
mathematical  properties  of  these  functions  have  not  been  fully  explored. 


CONCLUSIONS 

In  this  appendix  we  have  presented  the  error  behaviors  of  a binary 
DPSK  system  under  influences  of  intersymbol  interference  and  noise  correla- 
tions. The  graphically  presented  error  probability  curves  are  applicable 
to  system  performance  evaluations  when  SNR  imbalance  at  the  phase  detector 
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is  known.  Our  results  show  that  in  a DPSK  system,  the  error  probability 
does  not  depend  on  noise  correlation  when  the  symbol  probabilities 
are  equi-probable  regardless  whether  there  is  intersymbol  interference 


or  not.  The  noise  correlation  affects  the  error  probability  only  when 
the  message  symbol  probabilities  are  unequal. 


m 


J.  S.  LEE  ASSOCIATES,  INC. 


APPENDIX  I 

COMPUTER  PROGRAM  LISTING 
FOR 

ERROR  RATE  PERFORMANCE 
OF 

DPSK  OVER  A TERRESTRIAL  LINK 
WITH 

POWER  IMBALANCE  AND  CORRELATED  NOISE 
AT  THE  PHASE  DETECTOR 


1-1 


J.  S.  LEE  ASSOCIATES,  INC 


FORTRAN  IV  G LEVEL  il  NR  1 N DATE  » T90N0  08/53/5B  PAGt  0001 


0001 

IMPLICIT  REAL  *81  6-M.l  .0-2  1 

0002 

OlHtNSlGH  HS0USI16) 

00 03 

1000 

RC 60(5, 1, END-90001  RnO.lSUOB.NRSQ 

0004 

1 

fORMATWF 5.1 .131 

0005 

RE 60(5,2001) 18 SOL  I SIIMPJ  ,IAP*1 .NRSQI 

0006 

2001 

FORMAT | 16F5.1 ) 

0007 

WRI  TE (6#2>RMO.LS0DB 

0008 

2 

FORMAT!///'  RHO-  • ,F5.  1 ,5X,  'lAMB0A**2*'  ,F5.  1,  • OB'/ 

• • R**2  IOBI*»T23,#PICI*tT60*  *P  l£|R6RM  • ,T  60,  *P  1 E 1 SPACE  • • 1 

0009 

IS0»101**IIS0D»/ID1> 

0010 

L6HBDA-OSORTILSOI 

0011 

0PRH0*1  oo»rho 

001? 

OHRHO- 100-6H0 

0013 

OPPUH*UPAHO/0«RHO 

0014 

onuop*omrmo/oprmo 

0015 

OPf  InDO^L  6M806 

0016 

OPL2*OPL  *OPl  * 

0017 

OHL-1 .00-16HB06 

0016 

OML2-OML*OML 

0019 

00  900  f RS  * 1 , NRS  Q 

0020 

RSOOB-RSOLISI IRS  1 

0021 

R SO* ID 1 RS006/ lOll 

0022 

XARG»-0.25C0*RSQ*QPL2 

0023 

TP».?5DO*OMCCP*OML2*RSQ 

0024 

TN*  • ?5DO*QPOQM*OML  2*  RSQ 

0025 

LP*-.25DO*RSO*OPOOM#OPl2 

0026 

LN*-.25DO*RSO*OMOOP*OPL2 

0027 

FP  — . 5D0*R  SQ*OMl 2/OPRmO 

0028 

FN  — .5DO*RSQ*OML2/CMRmO 

0029 

XX-OE  API XARGI 

0030 

C6LL  DOSUH(TP«LP .fP.SP) 

0031 

PEP-.5DO*XX*OPRMO*SP 

0032 

CALL  DOSUM(TN,LN*f N,$N> 

0033 

PE  N*  • 5 DO* X X*GMRHQ*  SN 

0034 

PE*.5D0*<  PEP*PEN| 

0035 

WRITE (6,31 RSOCfc, PE .PEP .PEN 

0036 

3 

FORMAT (3X.F4 .1  ,T20,1P010.3,T40, 1P010.3.T60.  IPO  10.31 

0037 

900 

CONTINUE 

0038 

GOTO  1000 

0039 

9000 

STOP 

0040 

END 

FORTRAN  IV  G LEVEL  21  DOSUN  DATE  * 790*0  08/53/58  PAGE  0001 


0001 

SUBROUTINE  D0SU*MT,6RGl,ARGF'SUN| 

0002 

IMPLICIT  REAL *81 A-H,  1,0-21 

C H*Q 

TERM 

0003 

SUM* IDO 

0004 

C0EF-1D0 

0005 

TERM*  IDO 

0006 

MS  TOP*  501 •DABSf  T 1 

0007 

IF  lMSTOP.LT. 10?)HSTOP*J02 

0000 

DU  100  M-1,MST0P 

0009 

FM*M 

0010 

FMP'  -F  M*  1 DO 

0011 

COE  E *CO£F*(T/fM| 

001? 

0TERM*T£RN 

0013 

LGR-LAGERR* M, ARGLI 

0014 

OEO*ONEF  l(FM,FMPl,  ARGF1 

0015 

T£RM*C06F*LGR»0F0 

0016 

sum*sum#te  pn 

0017 

IF  I 0A6SI  TERM).GT.046S(QTERM) 1 GOTO 100 

0018 

IF  1 DABS  (TERM)  .IT  .06BS  ( SUM  1 *10- 6)  GOT  0200 

0019 

100 

CONTINUE 

00  20 

WRITE (6, 1 IMSTOP 

0021 

1 

FUR  MAT  1 • THE  SUM  DIO  NCT  CONVERGE  • MSTOP< 

0022 

STOP 

0023 

200 

RETURN 

0024 

END 
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FORTRAN  IV  G LEVEL  21  LAGERR 


0001 

REAL  FUNCTION  L AGERR  *81 n » X ) 

0002 

IMPLICIT  R E AL*6t  A— 0-2 1 

000  3 

IMNU.1,2 

0004 

1 

LAGERR* 100 

0005 

RE  TURN 

0006 

2 

S-100 

0007 

IFIU3.4,  3 

0006 

3 

T*  1 DO 

0009 

DO  100  H-1«N 

0010 

A»H 

0011 

T-T* ( n-H*1 DO  »/ ( A * A 1 

0012 

T*  T*  X 

0013 

T--T 

0014 

S-S*T 

0015 

100 

continue 

0016 

4 

L AGERR ■ S 

0017 

return 

0016 

END 

CA1E  * 79040  06/53/56 
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FORTRAN  IV  G LEVEL  21 


CNEF1 


0001 

REAL  FUNCTION  ONEF  1*  81  A,  B,  2 > 

0002 

IMPLICIT  R E AL*8 I A— H.  0-21 

0003 

T-  100 

0004 

5*  100 

0005 

AN1-A-1D0 

0006 

BMl-B-100 

0007 

00  100  !• 1 , 100 

0006 

C-l 

0009 

T-TA2/C 

0010 

T-T»IIAM1*CI/(BN1«CI) 

0011 

S*S»T 

0012 

IF  I0ABSI  T/SI  .LT.  10-81  GOT0200 

0013 

100 

CONTINUE 

0014 

WRITE (6 .1 ) A.B.2 

0015 

1 

format  i • if  11  • .o  15.8.  • . • ,015. 1 

(8  ACCURACY') 

0016 

2 00 

ONEFl-s 

0017 

return 

0016 

END 

DATE  ■ 79040  08/53/58 
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